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Other  systems.  Study  of  the  master  kinetic  equations  demonstra.ed  that 
approximate  analytical  solutions  exist  under  certain  conditions.  It  is  shown 
that  quantitative  prediction  of  the  vibrational  population  distribution  can  be 
determined  from  knowledge  of  just  a few  key  parameters  in  many  cases  of 
practical  interest.  *rhe  experimental  measurement  of  small  signal  gains  in  highl^ 
expanded  C0-N2  mixtures  confirmed  the  anticipated  large  gains  achievable  at  very- 
low  temperatures.  However,  unexpected  rotational  nonequilibrium  was  observed 
which  can  limit  the  power  available  from  P branch  laser  lines.  A summary  of 
efforts  to  provide  a simplified  method  of  computing  vibrational  excitation  rates 
in  a discharge  is  presented.  A modified  E/N  parameter  was  established  which  is 
shown  to  correlate  results  determined  in  a variety  of  gas  mixtures.  An 
experimental  study  of  a double  discharge  excitation  scheme  in  a supersonic 
C0-N2  flow  disclosed  that  the  secondary  pin  discharge  is  not  capable  of 
providing  a stable  discharge  with  independently  controlled  E/N  under  the 
conditions  tested.—  Considerable  effort  was  made  to  measure  the  characteristics 
of  the  plasma  obta'ihed  and  to  characterize  the  arcing  phenoma  observed.  An 
outgrowth  of  the  experimental  study  was  an  analysis  of  the  double  discharge 
ionization  scheme  with  particular  emphasis  on  those  characteristics  which  lead 
to  arcing  such  as  that  observed.  A small  CO  waveguide  laser  was  constructed  and 
tested.  Although  lasing  was  observed,  performance  was  far  below  that 
estimated  in  a theoretical  study  (discussed  in  an  earlier  annual  report).  An 
analysis  of  the  electrodes  after  several  hours  of  operation  and  other  operating 
conditions  indicated  that  iron  impurities  in  the  CO  probably  was  the  major 
cause  of  poor  performance.  Small  signal  gain  measurements  showed  unanticipated 
nonuniformities  in  the  lasing  medium. 
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I.  INTRODUCTION 


This  is  the  final  report  on  work  for  AFOSR  contract  AF  44620-73-C-0059 
entitled  "Anharmonic  Oscillator  Lasers".  Emphasis  of  the  work  perfo>^med  under 
this  contract  has  been  placed  on  fundamental  theoretical  and  experimental 
studies  of  various  phenomena  associated  with  diatomic  molecular  gas  flows  in 
which  anharmonici ty  plays  an  important  role  in  producing  or  maintaining  an 
inversion  which  is  suitable  for  lasing.  Most  of  the  work  concerned  carbon 
monoxide  gas  dynamic  or  electrically  excited  flows  because  these  systems  were 
capable  of  being  studied  in  our  laboratories,  had  considerable  potential 
practical  importance,  and  displayed  many  characteristics  common  to  other 
diatomic  gas  lasers.  In  several  cases  the  results  may  be  generalized  to 
other  systems. 

While  this  is  a final  report,  emphasis  is  placed  on  results  obtained 
during  the  past  year.  In  many  cases  the  reader  will  be  referred  to  earlier 
annual  reports,  published  literature,  and  theses  for  additional  discussion  and 
details.  The  report  is  organized  into  a theoretical  section  (II)  and  an 
experimental  section  (III);  however,  each  of  the  three  broad  areas  studied 
involved  bothe  experimental  and  theoretical  portions. 

These  three  areas  are: 

1.  Calculation  and  Measurement  of  Vibrational  Population  Distri- 
butions Under  Quasi -Steady  Conditions 

Study  of  the  master  kinetic  equations  (Section  II-B)  demonstrated 
that  approximate  analytical  solutions  exist  under  certain  conditions. 
It  is  shown  that  quantitative  prediction  of  the  vibrational  population 
distribution  can  be  determined  from  knowledge  of  just  a few  key  para- 
meters in  many  cases  of  practical  interest.  The  experimental  measure- 
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merit  of  small  signal  gains  (Section  III-B)  in  highly  expanded  CO-N^ 
mixtures  confirmed  the  anticipated  large  gains  achievable  at  very  low 
temperatures.  However,  unexpected  rotational  nonequilibrium  was 
observed  which  can  limit  the  power  available  from  P branch  laser 
lines. 

2.  Electron  Excitation  and  a Double  Discharge  Excitation  Scheme 

Various  aspects  of  this  important  problem  area  were  studied.  Sec- 
tion II. A presents  a summary  of  efforts  to  provide  a simplified  method 
of  computing  vibrational  excitation  rates  in  a discharge.  A modified 
E/N  parameter  was  established  which  is  shown  to  correlate  results 
determined  ina  variety  of  gas  mixtures.  An  experimental  study  of 

a double  discharge  excitation  scheme  in  a supersonic  CO-N2  flow 
(Section  III. A)  disclosed  that  the  secondary  pin  discharge  is  not 
capable  of  providing  a stable  discharge  with  independently  controlled 
E/N  under  the  conditions  tested.  Considerable  effort  was  made  to 
measure  the  characteristics  of  the  plasma  obtained  and  to  character- 
ize the  arcing  phenomena  observed.  An  outgrowth  of  the  experimental 
study  was  an  analysis  of  the  double  discharge  ionization  scheme 
with  particular  emphasis  on  those  characteristics  which  lead  to  arcing 
such  as  that  observed.  The  preliminary  results  of  this  study  are 
presented  in  Section  II. C. 

3.  CO  Waveguide  Laser  and  Time  Dependent  Electron  Pumping 

A Small  CO  waveguide  laser  was  constructed  and  tested  (Section 
III.C).  Although  lasing  was  observed,  performance  was  far  below 
that  estimated  in  a theoretical  study  (discussed  in  an  earlier 
annual  report).  An  analysis  of  the  electrodes  after  several  hours 
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of  operation  and  other  operating  conditions  indicated  that  iron  im- 
purities in  the  CO  probably  was  the  major  cause  of  poor  performance. 
Small  signal  gain  measurements  showed  unanticipated  nonuniformities 
in  the  lasing  medium. 
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II.  THEORETICAL  STUDIES 


II. A.  Vibrational  Pumping  by  Electron-Molecule  Collisions 

Knowledge  of  the  electron  distribution  function  is  critical  to  the 
determination  of  vibrational  pumping  rates  in  any  electrically  pumped  laser. 
When  the  vibrational  distribution  is  computed  in  the  standard  way,  by  inte- 
grating the  time  dependent  master  equations,  the  distribution  function  enters 
directly  into  the  rates  as  follows: 

30 

y*S^.j(c)f^(c)cdc 

0 

Here  N^.  is  the  number  density  of  the  i^^  vibrational  state  of  the  molecule 
Og  the  electron  number  density,  the  cross  section  of  the  i to  j 

transition  by  electron  impact  at  energy,  e and  fg(e)  is  the  electron  dis- 
tribution func'  For  the  quasi-steady  calculation,  the  determination  of 

the  specific  depends  on  the  total  quantum  pumping  rate  into  the 

vibrational  itaies.  This  net  pumping  rate  is  similarly  related  to  the 
electron  energy  distribution  function  when  the  pumping  mechanism  is  by 
electron  impact.  The  external  pumping  rate  representing  the  net  input  of 
quanta  into  vibration  is: 


OC 


where  = n„N  and  5,  = N./N  and  N is  the  total  number  density.  Since 

the  electron  energy  distribution  generally  varies  rapidly  in  the  energy 

range  of  interest,  the  approximation  of  the  distribution  function  by  a 

1 2 3 

Boltzmann  one  at  some  electron  temperature  is  inaccurate.  ’ ’ 
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The  most  commonly  used  procedure  in  laser  applications  is  to  com- 
pute numerically  the  electron  distribution  function  for  each  mixture  and 
E/N  of  interest.  Generally  only  collisions  with  molecules  in  the  ground 
vibrational  state  are  considered  in  the  electron  distribution  function  cal- 
culation even  when  superelastic  collisions  are  included  in  the  vibration 
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kinetics.  ’ The  objective  of  this  research  effort,  started  three  years 
ago,  has  been  to  a)  ascertain  the  effect  of  superelastic  collisions  on  the 
electron  distribution  function  and  thus  the  rates  and  b)  to  explore  al- 
ternatives to  the  conventional  strictly  numerical  computations. 

The  initial  formulation  including  the  superelastic  collisions, 
some  approximate  solutions,  and  initial  numerical  results  were  presented  in 
the  Annual  Report  1974-1975.^  The  stress  at  that  time  was  on  the  inclusion 
of  superelastic  collisions.  The  main  conclusion  is  that  below  about  1500°K 
effective  low-vibrational -level  vibrational  temperature,  superelastic  col- 
lisions play  virtually  no  role  in  determining  the  distribution  function. 

At  higher  effective  vibrational  temperatures  the  distribution  function 
slope  is  changed  in  the  vibrational ly  active  region  and  the  hign  energy 
"tail"  is  raised  as  reported  previously.  As  expected, the  fraction  of 
power  going  into  vibration  at  a fixed  E/N  decreases  with  increasing  vibra- 
tional temperature.  Because  even  in  highly  pumped  situations  the  zero  to 
one  vibrational  temperature  is  seldom  very  high,  the  exclusion  of  super- 
elastic effects  on  the  distribution  function  are  therefore  generally  not 
very  important. 

Subsequent  effort  was  devoted  to  the  development  of  a scaling 
law  that  would  eliminate  the  necessity  of  lengthy  numerical  computations 
for  each  specific  mixture  and  E/N.  A modified  E/N  parameter  arising 


naturally  from  the  analysis  was  found  to  scale  the  results  of  He,  CO, 
mixtures  over  a wide  rangt  of  mixture  ratios  as  long  as  electronic  ex- 
citation was  not  very  important.  Some  initial  results  were  presented  at 
the  28th  Gaseous  Electronic  Conference  in  Rolla,  Missouri,  in  October  1975. 
A more  complete  report  was  published  in  the  Proceedings  of  the  Princeton 
University  Conference  on  Partially  Ionized  Uranium  Plasmas  (NASA,  Sept. 
1976).  A preprint  of  the  latter  publication  wa^  attached  as  appendix  A 
to  the  Annual  Report  1975-76.^ 

The  main  conclusions  of  that  work  follow.  A modified  E/N  scaling 
parameter  which  can  be  interpreted  as  the  ratio  of  energy  gained  from  the 
electric  field  between  typical  inelastic  collisions  to  a characteristic 
energy  exchanged  in  such  an  inelastic  collision  arises  naturally  and  ex- 
plicitly from  the  analysis.  The  correlation  is,  however,  essentially  a 
local  one  in  energy  space,  i.e.,  it  is  the  local  logarithmic  slope  of  the 
distribution  function  that  scales  with  the  modified  E/N  parameter  based  on 
the  dominant  inelastic  process  in  that  energy  range.  For  the  cases  where 
vibrational  excitation  is  dominant  the  distribution  function  and  its  rele- 
vant moments  were  shown  to  correlate  to  within  about  10?^  for  a wide  range 
of  mixtures  of  CO,  N2  and  He.^  For  sufficiently  high  E/N,  electronic  ex- 
citation contributes  significantly  to  energy  exchange,  and  it  is  clear  that 
the  distribution  function  at  higher  energies  does  not  correlate  well  with 
respect  to  a parameter  based  solely  on  vibrational  collision  processes. 
Local  correlation  on  a parameter  based  on  electronic  cross-sections  is 
possible,  but  clearly  the  value  of  the  correlation  decreases  as  the  number 
of  correlation  parametersincreases.  However,  vibrational  excitation  rates 
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(the  important  results  for  laser  applications)  are  still  very  well 
correlated  by 


a 


■k  — 


E 1 


where  is  the  vibrational  energy  spacing,  & the  half  width  in  energy 

space  of  the  product  of  total  vibrational  cross  section  Q . 
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1 1 . B . Vibrational  Kinetics  of  CO  Laser  Systems 


Given  the  relevant  kinetic  rates  the  complete  master  kinetic 
equation  for  the  CO  vibrational  population  distribution  can  be  written 
down  and  solved  on  a digital  computer.  Such  a computer  code  was  written 

g 

at  Princeton  by  J.  Boulnois,  applicable  to  electric-discharge  excited  CO 
laser  systems.  In  addition  to  making  performance  calculations  for  the 
CO  wave-guide  laser,  this  code  was  used  to  generate  "exact"  quasi- 
equilibrium  solutions  for  comparison  to  other  approximate  solutions  ob- 

9 

tamed  using  analytical  techniques. 

The  rationale  for  pursuing  analytical  solutions  for  quasi - 
equilibrium  CO  vibrational  population  distributions  are  as  follows.  Even 
though  exact  solutions  can  readily  be  generated  by  a computer  code,  such 
numerical  solutions  inherently  require  explicit  specifications  of  the  rele- 
vant kinetic  rates.  At  the  present  time,  reliable  VV  and  VT  kinetic  rates 
are  available  only  at  low  quantum  numbers, and  the  required  qualitative 
features  of  the  exact  solutions  are  known  to  be  insensitive  to  the  assumed 
rates.  In  general,  the  first  few  levels  follow  approximately  a Treanor 
distribution,  the  very  high  levels  tend  toward  a Boltzmann  distribution, 
and  in  the  intermediate  levels  a relatively  flat  Plateau  distribution  pre- 
vails. It  is  of  interest  to  understand  the  physical  mechanisms  underlying 
these  observed  qualitative  behaviors  and  to  relate  quantiatively  the  popu- 
lation distributions  to  the  kinetic  rates. 

1 2 

Using  a differential  approximation  first  introduced  by  Brau, 


the  master  kinetic  equation  can  be  simplified  and  solved  analytically,  ex- 
pressing the  quasi-equi 1 ibri urn  CO  population  distribution  n^  explicitly  in 
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terms  of  the  kinetic  rates.  This  work^  was  presented  in  the  5^*^  Conference 

of  Chemical  and  Molecular  Lasers  (St.  Louis,  Mo.,  April,  1977)  and  has  been 

accepted  for  publication.  The  manuscriot  is  included  in  this  report  as  i 

1 

Appendix  B.  The  major  results  of  this  work  are:  j 

(i)  The  approximate  analytical  solutions  are  simple  and  are  sur- 
prisingly accurate  when  compared  to  "exact"  numerical  solutions  obtained 
by  the  computer  code. 

(ii)  At  given  T,  the  locus  of  the  boundary  between  the  Treanor 
and  the  Plateau  region  is  independent  of  the  kinetic  rates  assumed,  but 
is  only  a function  of  T,  a(anharmonicity)  and  9^^ j^(characteristic  vibra- 
tional temperature  of  CO).  An  explicit  formula  for  this  locus  NT  (T,i,9^^.j^) 
is  obtained. 


(iii)  At  given  T,  the  locus  of  the  PI ateau-Bol tzmann  knee  is  found 
to  depend  on  a certain  ratio  of  VT  to  VV  rates.  An  explicit  formula 
NT*  (VT,VV  rates)  is  obtained. 

(iv)  An  explicit  analytical  formula  for  the  Plateau  distribution 
N”^  is  derived  which  contains,  in  addition  to  N**  described  above,  an  inte- 
gration constant  which  represents  the  nominal  magnitude  of  the  Plateau. 

(v)  An  explicit  analytical  formula  is  obtained  relating  with 
^total  ’ total  vibrational  quanta  input  by  external  (electric 

discharge,  etc.)  excitation.  Thus,  once  is  known,  the  Plateau  dis- 

tribution can  be  simply  calculated.  Conversely,  given  the  Plateau  distribution 
(e.d.,  experimental ly( , the  value  of  can  be  deduced. 

(vi)  The  rate  of  vibrational  quanta  deactivation  by  VT  collisions, 
is  obtained  explicitly,  identifying  separately  the  contributions 


from  local  low  level  VT  collisions  and  high  level  VT  collisions. 


It  is 
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shown  thav.  under  weakly  pumped  conditions,  the  local  low  level  VT  colli- 
sions are  capable  of  deactivating  the  input  quanta,  and  the  resulting  CO 
population  distribution  will  not  exhibit  a Plateau.  Under  strongly  pumped 
conditions,  the  "excess"  input  quanta  which  can  not  adequately  be  locally 
deactivated  are  transported  to  high  levels  by  nearly  resonant  VV  collisions 
to  be  deactivated  by  VT  collisions  there.  Under  most  CO  laser  operating 

conditions,  the  local  low  level  contribution  is  small,  giving  rise  to  the 

VT 

interesting  result  that  D is  independent  of  VT  rates  but  is  dependent 
only  on  the  nearly  resonant  VV  rates. 

(vii)  The  anaiyiical  solution  shows  that  the  details  of  the  VV 
adiabaticity  factor  in  the  VV  rates  is  not  important  so  long  as  the  off- 
resonant  rates  are  small.  Only  the  second  moment  of  the  VV  adiabaticity 
factor  enters  in  the  approximate  solution  (via  NT*). 

(viii)  The  function  N|*  is  sensitive  to  the  magnitude  and  quantum 
number  dependence  of  the  VT  adiabaticity  factor. 

Since  the  quasi-equilibrium  CO  population  distribution  can  be  ex- 
plicitly related  to  the  VV  and  VT  kinetic  rates  via  the  analytical  theory, 
it  becomes  possible  to  infer  information  on  these  rates  particularly  at 

high  quantum  numbers  from  experimentally  measured  population  distributions. 

1 3 

In  collaboration  with  A.  Lightman  and  E.  R.  Fisher  of  Wayne  State  Univer- 
sity and  H.  Rabitz  of  Princeton,  we  have  successfully  inferred  some 
quantitative  information  on  CO-He  VT  rates  near  i = 25  at  150°K.  A draft 
of  this  work  is  included  with  this  report  as  Appendix  C. 


Before  discussing  this  work,  it  is  useful  to  review  briefly  a 
common  estimation  procedure  for  CO-He  VT  rates  at  high  quantum  numbers. 


n 
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Motivated  by  the  well  known  SSH  theory,  we  can  write 

= k'''''(04)  il^W(T)  (IIB.l) 

where  k'^'*'(i4)  is  the  i+1  ->■  i VT  rate,  is  proportional  to  the  square 

VT, 

of  the  transition  matrix  element,  and  (T)  is  the  VT  adiabaticity  factor. 
VT 

The  values  of  k (Of)  can  be  obtained  from  the  relatively  abundant  experi- 
mental data  {Z  '^'^=  1,  1). 

15  11 

Following  the  earlier  wcrk  of  Keck  and  Carrier,  Rich  et  al 
suggested  the  following  forms  for  and  aI^^: 


1+i 


A'f(T) 


(IIB.2) 

(IIB.3a) 


where 


3 ,,  1 

2 ^ ■ 3 ® 


-1^ 


3^ 


(IIB.3b) 


is  a curve-fit  formula  recommended  by  Keck  and  Carrier  as  a convenient  repre- 
VT 

sentation  of  A^.  (T)  computed  using  a simple  exponential  repulsive  interaction 
potential  via  the  SSH  theory.  The  parameter  is  given  by 


y.(T)  = (1-2  a.)  9 9'^'^/5T.  (IIB.3c) 

where  0'^^(°K'^)  is  a constant  proportional  to  the  reduced  mass  and  the  square 

of  the  interaction  length  L of  the  VT  collision.  In  this  rate  estimation 

VT  1 1 

procedure,  9 is  the  only  adjustable  parameter.  For  CO-He,  Rich  et  al 

VT  -1  ° 

suggested  9 = 0.114°K  , corresponding  to  L = 0.2  A.  For  large  y^ 
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(low  temperatures),  Ay"^(T)  simplifies  and  becomes  approximately 

= exp(3.)  (Rich  et.  al ) (IIB.4) 

where  

B = 4 ®vib 

At  low  temperatures  (T  < 300°K),  the  high  quantum  number  rates  estimated 
by  this  procedure  increase  with  i rather  steeply. 

Recently,  Verier  and  Rabitz^^  performed  extensive  numerical  cal- 
cui.'^ticns  for  CO-He  VT  rates  using  three-dimensional  quantum  mechanical 
close-coupling  methods.  Their  computed  1 -*■  0 rates  at  low  temperatures  agree 
within  an  acceptable  factor  with  available  experimental  data,  but  their 
high  quantum  number,  low  temperature  rates  disagree  substantially  with 
those  estimated  by  the  Rich  procedure.  As  reported  previously we  have 
correlated  the  Verter-Rabitz  rates  by  Eq.  (IIB.l),  obtaining  an  equivalent 
adiabaticity  factor  A^'^'ly^  ,T)  by  curve-fitting.  At  high  temperatures,  it 
is  found  that  Eqs.  (IIB.3)  correlated  the  computed  rates  adequately  if 
9'^^  is  chosen  to  be  0.271  °K”\  corresponding  to  L = 0.3  A which  is  in 
good  agreement  to  the  interaction  length  near  the  strongly  repulsive  part 
of  the  interaction  potential  used  in  the  calculations.  At  low  temperatures, 
a distinctly  different  correlation  results  and  is  given  below: 

[0.525  (1+4^)3 ~]  ij  (IIB.6a) 


(T  < 300  °K  only) 
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where 


1 


1 * (?|a)‘ 


exp 


(II.B.6b) 


In  Eq.  (IIB.6a),  the  exponential  factor  exp(-a6^^.j^/T)i  was  in- 
cluded because  of  detailed  balance  considerations.  The  theoretical  rates 
could  have  been  correlated  without  this  factor  (thus  yielding  a different 
set  of  curve-fitting  coefficients),  but  we  have  chosen  to  include  it 
because  it  formally  should  be  present  in  the  SSH  theory.  This  simple  cor- 
relation allows  us  to  succinctly  identify  the  distinctive  features  of  the 
Verter-Rabitz  rates.  These  are  (1)  the  presence  of  the  exponential  factor 
exp(-ct9^^^/T)i  in  Eq.  (IIB.6a)  which  renders  A^.  a much  more  gradual  i 

dependence  at  high  quantum  numbers*  and  (2)  the  presence  of  the  coeffi- 

VT 

cient  a^.  which  renders  A^.  a rather  raoid  i dependence  at  low  quantum 
numbers.  Of  particular  interest  is  the  unexpected  low  quantum  number 
behavior  of  a^  which  has  so  far  not  been  satisfactorily  interpreted  theo- 
retically. In  the  absence  of  experimental  data,  the  above  qualitative 
features  remain  to  be  substantiated  and  the  substantial  quantitative  dis- 
agreement with  the  Rich  estimated  rates  remairsto  be  resolved. 

From  a collection  of  experimentally  measured  CO  population  distri- 

1 3 

butions  obtained  by  Lightman  and  Fisher  in  a CO-He  electric  discharge 
laser  system,  we  have  inferred  CO-He  VT  rates  for  i = 25  at  150°K.  To 
our  knowledge,  these  are  the  first  experimental  VT  rates  at  such  high 
quantum  numbers. 

Our  results  as  given  in  Appendix  B indicate  that  (1)  the  in- 

*Note  that  S = 0.114)  and  0.525i(e'^'^  = 0.271)  are  reasonably  close 
numerically. 


ferred  rates  for  i = 25  and  150  °K  is  larger  than  the  theoretical  Verter- 
Rabitz  rates  by  a factor  of  about  7-8,  and  (2)  the  high  quantum  number  i 
dependence  of  the  inferred  is  gradual  and  is  consistent  with  the 
theoretical  Verter-Rabitz  predictions.  Since  the  CO-He  1^  0 VT  rates  are 
wel 1 -documented , the  observed  large  values  of  the  high  quantum  number  rates 
(with  gradual  i dependence)  suggests  that  at  low  quantum  numbers  A^  (or  a^. ; 
must  increase  rapidly.  Hence,  even  though  the  high  quantum  number,  T = 150  °K 
Verter-Rabitz  rates  substantially  underestimates  the  inferred  rates  quan- 
titatively, the  two  distinctive  qual i tati ve  features  of  the  Verter-Rabitz 
rates  are  consistent  with  the  inferred  rates.  Because  of  the  small  amount 
of  data  available,  the  accuracy  of  our  inferred  rates  is  rather  limited. 

But  the  data  reduction  technique  used  is  interesting  and  appears  to  be 
particularly  useful  at  high  quantum  numbers.  Additional  experimental  data 
are  clearly  needed  to  collaborate  and  refine  the  present  results. 


K 
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II. (C)  Further  Studies  on  Double-Discharge 

At  the  present  time,  no  self-consistent  and  complete  theoretical  model 
exists  for  the  double-discharge  configuration. 

The  governing  equations  for  a weakly  ionized,  continuum  plasma  flow 
are  the  following 

V • 4 = S 

V • = S 

• E = e(n.-ng) 

where  S is  the  net  volume  ionization  rate  and  and  are  number  current 
densities  given  by 

it  ■ "(U  - 
■ "ei  - 

and  is  the  velocity  of  the  neutral  flow.  and  are  diffusion  coeffic 
ients  and  and  are  mobility  coefficients.  Defining  j_  and  j_  by 

i = ig  ■ ii 

1 = + j i 

the  first  two  equations  in  (IIC.l)  can  be  rewritten  as 

V • i = 0 
7 • i = (1  + u)S 

and  Eqs.  (IIC.2)  become 

i = (ng-n.)u  - D.7(ng-n.) 

- (De-‘^i)'^"e  ■ ^"e'^e  ^ "i^i^i 


L. 


(IIC.l) 


(IIC.2) 


(lie. 3) 


(lie. 4) 


(lie. 5) 
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i = (^J%+n.  )u  - (D>uDg)Vn. 

^ (ni-ng)M4  + uDgV(n.-ng)  (IIC.6) 

If  the  local  Debye  length  of  the  plasma  is  small,  the  quasi-neutral  approxi- 
mation n^  - n^  = n can  be  invoked  , An  analysis  of  the  quasi -neutral  theory 
was  reported  in  the  last  Annual  Report.^ 

The  quasi-neutral  approximation  yields  a convection-diffusion  model 
for  the  plasma  density  n; 

u • Vn  = D,  7^n  + S - nV  • u (IIC.7) 

'—a 

where  0,  = (D.+uD„)/(l+u)  is  the  so-called  ambi-polar  diffusion  coefficient. 

a 16 

The  irrotational  electric  field  satisfies 

V • a E + (Dg-O.)V^n  = 0 (IIC.3) 

where  a = (u^+Ug)n  is  the  electrical  conductivity  (divided  by  e).  The 
coupling  between  Eqs.  (IIC.7)  and  (IIC.8)  occurs  through  S,  the  net  volume 
ionization  rate,  which  is  a sensitive  function  of  the  local  value  of  E/N 
where  N is  the  number  density  of  the  neutrals. 

The  quasi-neutral  formulation  accounts  for  convective  effects  and  is 
capable  of  describing  upstream  effects  (such  as  the  upstream  pin  discharge) 
on  the  developments  downstream.  However,  under  the  operational  conditions 
of  interest,  the  quasi-neutral  approximation  is  only  valid  in  some  limited 
local  regions  where  the  plasma  density  is  high;  it  is  clearly  not  applicable 
in  the  sheath  regions  adjacent  to  the  electrodes  and  in  the  region  beyond 
the  outer  edge  of  the  ambi-polar  diffusion  layer.  Because  of  its  non-uniform 
validity,  the  boundary  conditions  appropriate  for  the  model  are  not  yet  clear. 
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We  shall  confine  our  attention  to  the  cathode  sheath  and  the  sheath- 
diffusion  layer  interface.  Physically,  we  expect  the  normal  components  of 
and  to  be  continuous  at  the  interface  (which  shall  be  identified  by  the 
subscript  »).  The  quasi-neutral  model  yields: 


L = (i-iL 


/ r,  rs  \ 9n  V 1 

9n  1 

3yJac  ’ 


(D.+pDg) 


(I 

(I 


where  y is  the  unit  vector  normal  to  the  cathode  (and  therefore  the  thin  inter 
face)  and  E_  = -VV. 

In  the  cathode  sheath,  the  normal  component  of  £ (denoted  simply  by  E) 
is  expected  to  dominate.  Near  the  cathode  and  away  from  the  edge  of  the 
sheath,  both  convective  and  diffusive  effects  may  be  neglected.  Denoting  £ • 
by  j , Eq-  (lie. 5)  in  the  sheath  becomes 


j 


-(nebe 


n^.u^)E 


(I 


Denoting  ^ ^ ^q-  (IIC.6)  in  the  sheath  becomes,  using  the  Poisson's 

equation  in  Eqs.  (IIC.l): 


j = (n-ng)p.E  = 


^o'^i  r 3E 
e 3y  ■ 


Note  that  as  quas i -neutral ity  is  approached  near  the  sheath  edge,  Eq. 
becomes  untenable  while  Eq.  (IIC.ll)  can  still  be  justified. 


(I 

(lie. 12) 


Eqs.  (lie. 4)  in  the  sheath  can  be  written  as  follows 
j = constant 


dy 


(l+u)S  . 


(I 


ie.9) 

le.io) 


y 

le.ii) 

ie.i2) 


ie.i3) 


(lie. 14) 


Eqs.  (lie. 12),  (lie. 13),  and  (lie. 14)  govern  the  structure  of  the  sheath. 
We  shall  represent  the  net  volume  ionization  rate  by 


where  a(cm’b  is  the  first  Townsend  coefficient.  We  shall  represent  a by  the 
fol lowing  empirical  formula 


a = NB(N/E)  (lie. 15a) 

where  S(N/E)  is  to  be  determined  by  correlating  the  available  experimental 
data^^'^^  for  N2  which  has  been  collected  and  is  presented  in  Fig.  1.  It  is 
seen  that  for  low  N/E  the  data  can  adequately  be  represented  by 

e(N/E)  = a exp  (-b  I)  (IlC.lSb) 

where  a = 3 x 10"^^  cm^  and  b = 8.3  x 10'^^  volts-cm^.  For  high  values  of 
N/E  (N/E  -^0.7  X 10^^),  the  data  has  very  wide  scatter. 

We  can  rewrite  Eqs.  (IIC.12)  and  (IIC.14)  in  the  following  non-dimension- 
al i zed  form: 


E ^ = B J 
dn 


(u3)b 


where 

E = -E/bN,  3 = j/j  , 

n = aNy,  B = S(E) 

and  B is  the  non-dimensional  normal  current  density 
B = 


(lie. 16) 
(lie. 17) 


(lie. 18) 


ab^  £q  u^-  N’ 


(lie. 19) 


which  is  independent  of  n.  Dividing  Eqs.  (IIC.16)  by  (IIC.17)  and  integrat- 
ing once,  we  obtain 


3(§)dC  = B J - In  (1+j; 


(lie. 20) 


where  is  an  integration  constant  to  be  determined. 

The  above  mobility-dominated  theory  cannot  be  matched  directly  with  the 
convective-diffusive  quasi -neutral  theory  because  the  two  models  have  no 
region  of  overlapped  validity.  However,  it  appears  physically  reasonable  to 
require  that  the  sheath  solution  should  become  quasi -neutral  asymptotically 
as  n - This  ^ hoc  boundary  condition  imposes  = 0 and  J(rj==°)  = 0.  At 
the  surface  of  the  cathode,  we  shall  assume  that  the  electron  (secondary) 
emission  current  is  related  to  the  incoming  ion  current  by  the  second  Townsend 
coefficient  y:  . - 


Jp.  ■ "yJ-i 


J = J 


= blJ. 

cathode  y+1 


(lie. 21 ) 


at  n = 0-  Although  y can  be  as  small  as  0(10'^)  for  molecularly  clean  con- 
ductor surfaces,  for  common  engineering  surfaces  y = 0(10"^)  is  more  typical. 
Eqs.  (IIC.16)  and  (IIC.17)  have  been  integrated  numerically  using 

A 

"^cathode  ” given  by  Eq.  (IIC.15)  for  various  values  of  B. 

Figs.  2,  3,  and  4 show  the  distributions  of  E,  the  normalized  electron  and 
ion  densities,  and  the  normalized  voltage. 

As  mentioned  previously,  these  solutions  are  not  valid  near  the  edge  of 
the  sheath  away  from  the  cathode  because,  first  of  all,  the  data  for  3(E)  is 
not  reliable  at  low  values  of  E,  and  secondly  the  effects  of  convection 


(particularly  of  ions)  and  perhaps  of  diffusion  are  not  negligible.  However, 
Fig.  3 shows  that  the  sheath  solutions  approach  quasi-neutrality  rather 
rapidly,  and  the  nearly  neutral  normalized  plasma  density  varies  with  n very 
slowly.  Fig.  4 shows  that  in  the  nearly  neutral  region  the  normalized  vol- 
tage also  varies  slowly.  Fig.  5 shows  the  "sheath  voltage  drop"  between  the 

^ •'I 

cathode  and  a point  at  the  "edge"  of  the  sheath  with  J = 10  as  a function 
of  B.  It  is  seen  that  this  sheath  voltage  drop  varies  non-monotonical ly  with 

B . 

The  sheath  analysis  presented  here  is  intended  as  one  of  the  building 
blocks  toward  a self-consistent  and  complete  plasma  theory  for  the  double- 
discharge configuration.  The  sheath-diffusion  layer  matching  procedure  used 
is  an  ad  hoc  one  and  further  refinements  are  necessary.  The  region  beyond 
the  upper  edge  of  the  diffusion  layer  where  quasi-neutrality  also  fails  has 
not  yet  been  treated.  The  modeling  of  the  net  volume  ionization  rate  is 
rather  simple-minded  and  recombination  and  other  kinetic  processes  have  been 
ignored.  Much  further  work  remains  to  be  done  in  this  area. 


III.  EXPERIMENTAL  STUDIES 


III. A.  Double  Discharge,  Experiments  and  Modelling 

(1)  Introduction 

The  Princeton  electric  discharge  channel  is  a double  discharge 
scheme  for  the  production  of  stable,  distributed  plasmas  in  a supersonic 
(M  ^ 4)  stream  of  diatomic  nitrogen  and  nitrogen-carbon  monoxide  jnixtures. 
This  research  is  motivated  by  a desire  to  produce  electcically  pumped, 
near  infrared,  carbon  monoxide  lasing  action  from  stable,  large  volume 
plasmas  in  which  electron  number  density  and  average  energy  are  independent 
ly  and  externally  manipulated.  The  envisioned  devices  would  operate  con- 
tinuously and  be  of  high  output  power.  Our  research  efforts  are  aimed  at 

gaining  a physical  understanding  of  the  double  discharge  stabilization  tech 
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nique  introduced  by  Blom  and  Hanson,  while  attempting  to  determine  if 
present  devices  can  be  scaled  up  in  size  and  power  loading. 

Following  a brief  description  of  the  apparatus  (2),  summaries  will 
be  made  of  work  done,  results  obtained  and  hypotheses  entertained  in  each 
of  five  topics:  fluid  mechanics  studies  (3),  general  electrical  character- 
istics (4),  analysis  of  arcing  instability  (5),  preliminary  model  of  the 
plasma  stream  (6)  and  experimental  determination  of  plasma  properties  (7). 
This  report  will  concentrate  on  developments  since  the  completion  of  the 
Annual  Progress  Report  of  1975-76.^ 

(2)  Apparatus:  Dimensions,  Design  and  Description 

The  schematic  of  the  discharge  channel  (Figure  6)  plus  the  follow- 
ing tables  and  comments  should  adequately  describe  the  discharge  channel. 


TABLE  III  A.l 


Position 

Reference  Length 
(From  Throat) 

Height  Width 

cm  cm 

X-area 

cm2 

Throat 

0 

2.54  .635 

1.61 

Discharge  Volume 

Leading  Edge 

21 .90 

2.55  6.40 

16.31 

Trailing  Edge 

72.71 

2.93  6.40 

18.76 

TABLE  II] 

[ A. 2 

Nominal  Mach  Number 

3.93 

Area  Ratio 

10.08 

Type 

Contoured,  2D 

Length,  Throat  to 

19.69  cm 

Parallel  Section 

TABLE  III  A. 3 

ELECTRODES 

Item 

Geometry 

Materials 

Discharge  Plates 

(50.80  cm)  long,  span  full  width. 

Copper 

rounded  leading  & trailing  edges 

(latter  exposed) 

Pins 

.0635  cm  dia. 

Flush  with  cathode 

Tungsten 

surface  7 rows 

cf  5 across  the 

width 

Pin  Insulators 

Flush  with  cathode  surface 

Boron  Nitride 

( .2381  cm  dia. 

) 

The  first  row  of  pins  is  located  1.27  cm  from  the  leading  edge 
of  the  discharge  plates  with  four  successive  rows  at  2.54  cm  spacings.  Row 
five  and  six  are  spaced  at  7.62  cm  intervals.  The  center  pin  of  each  row 
is  on  the  tunnel  centerline,  the  rest  at  1.27  cm  intervals  on  either  side. 

(3)  Fluid  Mechanics 

It  is  useful  to  review  our  estimates  of  general  flow  properties 
computed  with  one  dimensional  isentropic  gasdynamic  theory.  Assuming  choked 
flow,  one  can  derive  the  following  expressions  for  mass  flow  in  the  dis- 
charge channel : 

TABLE  III  A. 4 
MASS  FLOW 


T 

0 

m (gm/sec) 

m( 1 bm/sec) 

0 

^o> 

K 

44.205  P^/VT~ 
o'  0 

.09740  P /vT~ 
0 0 

300° 

K 

2.552  P 

0 

5.627  X 10"^  P 

0 

where ; 

1) 

Pq  in  psia,  stagnation  pressure. 

2) 

(multiply  above  by  14.696  for  P^  in  ATMS 

.) 

3) 

T is  stagnation  temperature  in  degrees 

Kelvin 
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For  M = 3.93  the  number  density  and  pressure  can  be  calculated 
from  the  following  formulas: 


TABLE  III  A. 5 


TEST  SECTION  DENSITY 

T- 

N#/cm^ 

Cavity  Pressure 

Equivalent  Pressure 

0 

at  cavity 

at  density  N and 

(Pq  in  ATMS) 

temperature 

temperature  T = 300° K 

20 

2.155  X 10^^ 

Po/^o 

P = 22.474  P (^) 
0 ^0 

P = 5.497  P„ 
0 

300 

7.185  X 10^^ 

P 

0 

P = 22.474  P 

0 

where  P in  TORR 

P in  ATM. 
0 

Pitot  surveys  were  made  and  reported  in  the  Annual  Report  of  1975- 
76.^  These  data  indicated  that  while  there  is  a region  of  reasonably  uniform 
flow  a very  complicated  flow  field  exists  in  the  discharge  channel  which 
cannot  be  fully  analyzed  without  additional  sophisticated  flow  measurement 
techniques.  The  value  of  the  results  did  not  warrant  such  an  extensive 
and  expensive  investigation. 

The  boundary  layers  probed  have  been  thick,  ranging  from  .6  to 
1 cm  along  the  tunnel  sidewalls.  Boundary  layers  along  the  electrodes  are 
not  of  uniform  thickness  and  appear  to  bulge  near  the  center  line  of  the 


r 
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tunnel.  Pitot  profiles  exhibit  noticeable  peaks  and  valleys  which  make 
any  simple  identification  of  boundary  layer  edges  difficult.  It  is  also 
probable  that  oblique  wave  patterns  exist  in  both  planes  along  the  flow 
axis.  Our  nozzle  contour  could  be  inaccurately  machined  creating  an  ex- 
pansion fan;  or  the  break  from  the  contoured  section  to  parallel  walls 
might  be  too  sharp,  creating  shock  waves.  In  the  plane  perpendicular 
to  the  nozzle  expansion  we  might  have  a shock  pattern  induced  by  rapid 

electrode  boundary  layer  growth  or  leakage  across  the  throat-top  and  bottom 
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wall  seals;  these  effects  are  discussed  by  Simons"  and  Director. 

As  the  discharge  channel  is  a long  duct  of  small  cross-sectional 
area,  viscous  effects  will  be  important  even  at  supersonic  velocities.  The 
ratio  of  electrode  length  to  hydraulic  diameter  is  about  13.  An  estimate 
of  the  level  and  variation  of  flow  variables  can  be  obtained  on  the  basis 
of  a one-dimensional  frictional  flow  analysis.  Using  the  estimated  Mach 
number  at  the  leading  edge  of  the  electrode  of  = 3.82,  the  computed  flow 
preoerties  at  the  trailing  edge  are  presented  in  Table  III  A. 6 below. 

TABLE  III  A. 6. 

VISCOUS  EFFECTS 


M 

P/^L.E. 

°^"L.E. 


^/\.E. 


P /P  , r 
O'  oL.E. 


Leading  Edge  (L.E.) 
3.83 
1 
1 

1 

1 


Trailing  Edge  (T.E.) 
2.77 
1.72 
1.11 

1.55 

0.37 
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This  estimate  of  M(T.E.)  is  not  inconsistent  with  one  that  can 

be  deduced  from  pitot  and  wall  static  pressure  data  in  conjunction  wi th  this 

estimate  of  Pg(T.E.).  We  predict  about  an  in  increase  in  density  and  an 

4 4 

equivalent  drop  in  velocity  from  about  6.7  x 10  cm/sec  to  6 x 10  cm/sec. 

(4)  General  Electrical  Characteristics 

A general  circuit  diagram  of  the  double  discharge  is  provided  in  Fig- 
ure 7.  The  desired  operating  feature  in  the  double  discharge  design  is 
the  ability  to  independently  set  plate  voltage  and  current  at  a given 
(or  cavity  density),  that  this  has  been  accomplished  can  be  seen  in  Figure 
8 . It  should  be  noted  that  the  total  variability  of  (E/N)  (electric  field 
to  total  number  density)  for  this  device  is  quite  small,  ranging  within 
2 X 10"^^  v-cm“  to  4 X 10*^®  v-cm^. 

The  major  parameters  affecting  double  discharge  operation  are 
^o’  ^0  current),  cathode  current  density  and  arcing  rate.  The 

most  noticeable  effectson  discharge  operation  are  exhibited  by  changes  in 
the  parameters  P^  (stagnation  pressure  which  is  linearly  proportional  to 
cavity  pressu'^e,  density  and  mass  flow)  and  I^.  At  any  given  there  is 
a minimum  necessary  in  order  to  have  arc  free  operation  (or  specified 
arcing  frequency)  with  a desired  level  of  main  current,  I. 

As  Iq  increases,  the  amount  of  main  current  that  can  be  delivered 
without  arcing  increases;  this  effect  is  duplicated  if  P^^  decreases  and 
is  maintained. 

The  onset  of  arcing  is  gradual,  and  arcing  frequency  increases  for 

decreasing,  "insufficient"  I or  increasing  "excessive"  P . 

^ 0 0 


An  arc  formed 
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between  the  electrode  plates  must  be  swept  out  of  the  discharge  volume 
before  another  can  form.  Arcs  can  be  visually  detected  in  pinkish 

peachy  flames.  Frequent  arcing  manifests  itself  on  main  voltage  versus 
main  current  plots  as  a region  of  negative  slope. 

A measurement  of  arcing  rate  versus  average  anode  current  was  made 
at  one  and  fixed  pin  geometry  (Figure  9).  The  gradual  arcing  thresholds 
are  clearly  visible  as  is  the  effect  of  pin  current.  A sample  error  bar  is 
shown  on  our  values  of  average  current  in  the  linear  portions  of  the  graph. 
This  error  diminishes  rapidly  in  the  regions  of  low  arcing  rate.  Further 
discussion  of  arcing  is  given  in  the  following  sections. 

(5)  Arcing 

As  mentioned  in  the  previous  section,  the  major  instability  of  the 
double  discharge  is  the  formation  of  arcs.  Experiments  were  done  in  which 
it  was  shown  that  solitary  arcs  are  formed  and  are  swept  out  of  the  dis- 
charge volume  at  the  flow  speed  before  a recurrence.  Photomultiplier  tubes 
were  used  to  monitor  arcs  as  they  passed  one  location  in  the  tunnel  in  an 
experiment  described  in  the  Annual  Progress  Report  of  1975-76.^ 

In  an  effort  to  check  and  refine  the  results  of  the  photo-tube 
data  another  experiment  was  performed.  One  photo-multiplier  observed  two 
posistions  in  the  discharge  channel  separated  by  a known  distance.  An 
oscillogram  of  the  photo-tube  signal  would  show  two  peaks  representing  an 
arc  passing  each  of  the  monitored  positions.  This  direct  measurement  cf 

4 

arc  velocity  yields  a value  of  5.4  x 10  cm/sec. 

A schematic  of  an  oscillogram  is  given  in  Figure  10.  The  peaks 
in  the  oscillogram  are  of  different  magnitude  because  each  was  transmitted 
by  a different  leg  of  an  uncalibrated  optical  system.  The  RC  delay  of 
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the  photo-tube  circuitry  was  measured  and  minimized  sufficiently  that  our 

-4 

ultimate  detectability  of  "arc  thickness"  was  .013  us  or  7 x 10  cm 

4 

moving  at  the  flow  speed  of  5.4  x 10  cm/sec.  The  half-width  of  these 
signals  is  about  60  us,  indicating  a maximum  arc  thickness  of  3.2  cm.  The 
same  data  could  however  be  explained  by  a much  thinner  sheet  inclined  with 
respect  to  the  optical  axis. 

(6)  Plasma  Model,  A Synopsis 
1.  Governing  Equations 

The  Debye  length  likely  to  be  encountered  in  the  double  dis- 

-3 

charge  is  of  the  order  of  10  cm.,  hence  a continuum  theory  of  ambipolar 
diffusion,  convection  and  production  of  the  plasma  in  the  bulk  of  the  test 
section  should  be  appropriate.  Assuming  singly  charged  ions,  negligible 
attachment,  steady  state,  quasi -neutral i ty  and  small  ionization  fraction, 
governing  equations  for  the  plasma  flow  can  be  derived  as  in  Ref.  26. 

The  application  of  this  formalism  to  modelling  of  the  double-dis- 
charge configuration  was  described  in  the  last  annual  report^  and  some 
simplified  solutions  were  presented.  For  convenience  some  results  from  that 
report  for  the  infinite  space  and  semi -infinite  space  solutions  due  to  a 
single  source  term  are  shown  in  Figure  11a  and  11b  respecti vely . We  expect 
that  the  semi-infinite  space  solution  is  the  more  relevant  to  our  configuration 
and  is  indicative  of  the  primary  spatial  variation  of  number  density. 

The  analytic  expression  of  Equation  IV  (F.17)  in  Reference  7 can 
be  used  as  a guide  in  anticipating  plasma  density  distribution  data.  If  we 
consider  the  case  of 


x 
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where  (x-x^)  is  the  streamwise  distance  and  (y-y^)  and  (z-z^)  are  the 
transverse  distances  from  the  source,  we  can  find  an  approximate  expression 
for  the  variations  of  the  number  density  n with  z and  y with  a given  cross 
section; 


n « exp 


- 3 ((y-y„ 


where  k = u/D  ,the  flow  velocity  u divided  by  the  ambi polar  diffusion 

d 

coefficient  D .and  5 = cs/D  where  a is  the  ionization  rate. 

d d 

From  this  expression  we  would  anticipate  data  at  a given  cross 

2 

section  when  plotted  semi-logarithmically  versus  the  coordinate  (y-y^) 

2 

+ (z  - Zq)  to  be  linear  with  negative  slope,  and  when  plotted  versus 
y or  z to  be  parabolic. 

An  estimate  of  the  ambi polar  diffusion  coefficient  can  be  made 
from  the  slope  of  the  lines  just  described.  If  ^ » a as  appears  to 
be  true  at  least  in  certain  locations  in  the  discharge  volume  slopes  of 
ln(n)  versus 


2 2 

(y-y^)  + should  be  proportional  to 


- 1 Ji  = . u 

] x-x  1 4 4 1 x-x„ I D 

I o'  ' o'  a 

With  u known,  D can  be  estimated  from  the  data. 

d 

If  all  data  at  a given  cross  section  cannot  be  correlated  by  a 
single  straight  line,  then  the  model  may  still  be  valid  but  require  a dis- 
tributed source  instead  of  the  single  concentrated  one.  As  the  model  equa- 
tion is  linear,  one  can  experiment  with  source  distributions  in  order  to 


analytically  duplicate  experimental  results.  Electron  and  ion  distributions 
between  the  electrodes  at  x = are  not  known  because  of  low  plasma  density 
and  resulting  non  quasi -neutral i ty.  There  is,  therefore,  a certain  amount 
of  freedom  in  the  initial  conditions  which  can  account  for  some  discrepancy 
between  the  model  and  the  data.  If  the  resulting  values  of  are,  however, 
reasonable  and  the  model  predicts  distributions  at  several  cross  section 
locations  it  is  most  likely  that  the  major  assumptions  are  correct. 

A test  can  be  made  to  determine  if  volume  ionization  does  or  does 


not  play  a role  in  the  bulk  of  our  discharge  volume.  It  is  well  known  that 

the  dependence  of  ionization  rate  (a)  on  the  electric  field  E is  extremely 

2 

rapid,  being  much  stronger  than  that  of  D , hence  k /4  is  almost  independent 

d 

of  E compared  to  the  rapid  variation  of  a. 

One  can  measure  density  distributions  at  one  cross  section  of  the 
tunnel,  first  with  sustainer  plates  active  to  enhance  S and  then  with 
sustainer  plates  off,  at  one  value  of  pin  current.  Plots  of  ln(n)  versus 
(y-y^)^  + {z-Zq)  will  have  slopes  proportional  to  — a,  if  this  drop- 
off in  density  is  much  slower  in  the  active  sustainer  case  than  in  the 
passive  sustainer  case  and  the  level  of  density  is  higher  at  a given  location, 
then  volume  ionization  will  have  been  shown  to  be  a significant  effect  in 
producing  the  plasma  cloud. 


The  most  desirable  measurements  that  could  be  made  in  a double- 


discharge plasma  would  be  spatially  resolved  electron  number  density  and 
potential  measurements.  This  requires  the  design,  construction  and 
successful  employment  of  Langmuir  probes. 
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The  probing  scheme  developed  during  the  past  year  is  shown  in 

Figure  12.  As  our  electrons  have  an  energy  of  about  1 ev  an  electric 

probe  in  the  plasma  should  saturate  when  biased  by  a few  vclts  from  floating 

potential.  A measurement  of  this  saturation  current  can  be  related  to 

plasma  density.  Our  technique  is  to  send  a five  volt  200  Hz  sine  wave 

through  a d.c.  blocking  capacitor  to  the  probe.  Circuit  components  are 

sized  so  that  the  probe  to  plasma  voltage  drop  is  practically  the  entire 

five  volts.  The  current  at  200  Hz  is  monitored  on  the  low  voltage  side 

of  the  capacitor  by  measuring  the  differential  rms  voltage  across  a 

precision  resistor  with  a lock-in  amplifier.  The  frequency  of  200  Hz 

was  chosen  because  the  plasma  is  relatively  quiet  there.  This  was  checked 

with  a spectrum  analyzer.  In  the  range  of  100-300  Hz  we  see  only  white 

noise  at  about  1.6  volts;  our  lock-in  amplifier  has  no  trouble  monitoring 

our  200  Hz  signal.  These  readings  of  saturation  current  can  be  plotted 

versus  position  to  produce  a relative  plasma  density  map.  Figure  13  is 

a olot  of  such  signals  versus  — (.5  = tunnel  centerline)  at  a given 

^total 

y locations  from  the  cathode  surface.  Figure  14  is  a cross  plot  from 

Figure  13  which  shows  contours  of  constant  signal  level  at  a given  cross 

section.  The  data  spans  two  orders  of  magnitude  in  density. 

2 

An  estimate  of  D from  the  data  is  24.2  cm  /sec.  Using  the  fact 

cl 

that  electrons  are  much  hotter  than  ions  in  this  type  of  plasma  and  from 
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data  on  ionic  mobility  coefficients  we  can  estimate  D^.  Using  the  for- 

a 


mu  I a : 


kTe 


e 


with  — - = 1 volt, 

e 

E = 500  v/cm, 

p ? 73  °K  12.2  torr, 

Pequivalent  ^ 

2 

we  find  D,  = 26  to  50  cm  /sec  in  nitrogen  depending  on  whether  atomic  or 

d 

molecular  ions  are  involved.  Though  our  data  is  limited  we  see  that  the 
model  has  some  validity.  A reasonable  D,  also  appears  to  describe  che 

a 

amount  of  spanwise  spreading  observed  near  the  cathode. 

It  appears  that  some  source  distribution  instead  of  just  a delta 
function  at  x = (i.e.  upstream),  is  more  likely  to  be  consistent  with 

the  present  data. 

In  order  to  arrive  at  an  absolute  number  density  estimate  from 
our  measurements  of  saturation  currents  we  proceed  as  follows: 

I = enuA 

e = 1.602  20'^^  coul. 
n = number  density 
u = flow  velocity 

A = cross  section  of  a stream  tube  whose  bounding 
streamlines  traverse  the  probe  bow  shock  and 
graze  the  edge  of  the  probe  sheath. 

We  say  that  the  quasineutral  plasma  is  swept  downstream  at  the  flow  velocity 
flow  velocity  and  that  all  particles  in  the  stream  tube  are  carried  into  an 
electric  sheath  where  they  are  retrieved  by  the  probe.  Our  problem  is  in 
estimating  A.  Most  relevant  theories  of  plasma  probing  surveyedseem  to 


give  estimates  of  probe  sheath  thickness  within  a factor  of  2 or  3 of  each 
other  and  with  a slew  dependence  on  n.  For  our  experiment  we  expect  the 
ratio  of  sheath  thickness  to  probe  diameter  to  be  smaller  than  or  nearly 
equal  to  1.  Therefore  as  a first  estimate  we  say  that  A is  approximately 
equal  to  the  frontal  area  of  the  probe  wire  and  is  a constant  value  where- 
ever  we  make  our  measurements.  In  actuality  A increases  as  n decreases 
reducing  spatial  resolution  in  regions  of  lower  plasma  density  and  signal 
level . 

Between  the  value  of  noise  level  and  the  maximum  signal  our  in- 
strument can  assimilate  (the  range  in  Figure  13)  we  estimate  an  ability 
to  detect  densities  from  10^^  to  7 x 10^^  cm”^. 

As  a final  note  we  should  say  that  we  have  measured  a floating 
potential  rise  of  about  230  volts  within  0.1  cm  of  the  cathode  surface. 

(8)  Summary 

Brief  reviews  were  made  of  the  results  of  earlier  flow  and  ter- 
minal electrical  characteristic  measurements,  as  well  as  of  the  arc 
propagation  and  plasma  flow  models. 

Data  was  presented  on  the  dependence  of  arcing  rate  on  discharge 
parameters  and  the  arc  propagation  model  was  checked  with  a direct  measure- 
ment of  arc  speed  and  width.  Cathode  fall  voltages  have  been  measured  at 
some  locations  and  found  to  be  about  200  to  230  volts. 

A successful  plasma  probing  technique  was  developed  and  relative 
plasma  density  profiles  and  contours  have  been  mapped  out  at  the  middle 
window  port  with  one  pin  and  the  sustainer  plates  active.  Comparison  of 
data  with  theory  indicates  a need  to  improve  our  modelling  of  the  source 
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terms  but  does  not  negate  as  yet  the  foundations  of  the  theory.  More  data 
is  needed  to  further  test  the  soundness  of  our  theoretical  model  and  to 
elucidate  the  role  of  volume  ionization  in  the  discharge  volume. 
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1 1 1. 3.  CO-Gain  Measurements 


Small  signal  gain  measurements,  made  in  a low  temperature  CO-N2 
mixture  expanded  through  a hypersonic  nozzle,  were  reported  in  the  last 
annual  report.^  Static  temperatures  were  varied  over  a narrow  range  (30, 

35,  41  °K)  by  altering  nozzle  throat  diameter.  Stagnation  conditions, 

2000  psia  and  2000  °K,  and  gas  mixtures,  C0:‘J2  = 15:85,  were  maintained 
fixed  throughout. 

Attentio'n  here  is  focused  on  the  M = 18.5  data  (T.  = 30.3°K, 

' trans 

16  3 

H|^^_Co  = 1.2  X 10  cm"'^).  Specifically,  this  discussion  concerns  the 

detection  of  rotational  nonequilibrium  deduced  from  the  vibrational  popu- 
lations computed  from  measured  small  signal  gains.  Details  of  the  three 
pass  optical  system  •mployed  to  measure  gain  may  be  found  in  Refs.  7 and  28. 
The  measured  gains  are  presented  graphically  in  Ref.  7 and  are  tabulated  in 
Ref.  28.  Vibrational  populations  corresponding  to  the  values  of  gain 
measured  for  = 30.3  are  shown  in  Figure  15.  These  populations  were 

derived  from  the  small  signal  gain  equation  for  a Doppler  broadened  line 
using  the  equation 
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where  A^g  = spontaneous  transition  probability  = 30.3  sec 
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'10 

K 
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= 2143.2  cm' 
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= 1.38  X 10’’^  ergs/“K 
= 4.648  X 10"^^  grams. 
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From  Reference  29 

a 

^ ( (l-2x  V-x  ](l-2x  V+x J 1 

j (l-3Xg)(l-XgV)  i’ 

where  x^  = anharmonicity  factor  - .006126, 
and 

'^rot  " ^ (20+1)  exp  {- 

J”  1 


V,V-1 
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l-2x  V 
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(III.B.2) 


(III.B.3) 


where  9 = characteristic  rotational  temperature  = 2.78  °K.  The  symbol  T 
is  used  to  denote  the  rotational  temperature  as  defined  from  the  relative 
population  of  rotational  states  0-1  and  0,  in  the  most  general  case. 

The  vibrational  populations  shown  in  Figure  15  correspond  to  solu- 
tions of  Eq.  III.B.1  using  T = = 30  °K  and  pairs  of  gain  values  for 

two  P branch  transitions  on  the  same  vibrational  band.  The  particular  pair 
of  P branch  gains  used  to  compute'  are  denoted  on  the  figure.  The  open 
symbols  denote  values  of  n,y  computed  from  gains  on  the  V band,  whereas  the 
closed  symbols  were  computed  from  gains  on  the  V+1  band.  Three  features 
are  to  be  noted  in  Figure  15. 

1.  There  is  a surprising  amount  of  scatter  in  the  populations, 
indicating  that  the  populations  predicted  are  greatly  in- 
fluenced by  the  particular  pair  of  gain  values  used  in  the 
computation.  This  is  disturbing  since  according  to  the  theory 
and  associated  assumptions  these  points  should  lie  on  top  of 
one  another. 

2.  A general  tendency  exists  that  the  populations  ccirputed  from 
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lower  rotational  states  are  lower  than  those  computed  from 
higher  rotational  transitions. 

3.  The  V=6  populations  computed  from  V=7  band  data  are  approxi- 
mately a factor  of  four  larger  than  those  derived  from  V=6 
data. 

It  was  postulated  in  Reference  7 that  the  major  source  of  the  scatter 

arose  from  the  possibility  of  rotational  nonequilibrium.  Initially,  values  of 

T greater  than  were  employed  in  Equation  III.B.l,  with  the  result  that 

the  scatter  reduced  somewhat  as  T increased  until  T exceeded  35  °K 

(T.  = 30  °K),  at  which  point  the  scatter  increased  as  T was  raised  further. 

uranS 

This  was  not  only  unsatisfactory  as  far  as  collapsing  values  of  ny  to  a series 
of  single  points,  but  it  is  inconsistent,  we  believe,  to  assume  that  if  non- 
equilibrium exists  between  rotational  and  translational  modes,  the  rotational 
populations  can  be  described  by  a simple  Boltzmann  distribution.  (For  J ^6 
the  rotational  spacing  2J0  exceeds  the  translational  temperature  at  this  low 
temperature.)  The  scatter  in  Figure  15  is  greater  than  car  be  attributed 
to  errorsin  gain  measurement  particularly  for  the  V=7  and  8 bands.  Gains  were 
derived  from  an  average  of  three  separate  measurements,  which  for  all  but  the 
smallest  gains  varied  by  less  than  + 10%. 

There  is  experimental  evidence  of  another  sort  that  indicates  rota- 
tional nonequilibrium  exists  under  conditions  similar  to  those  studied  here. 

In  the  1960's  several  attempts  to  determine  gas  temperature  using  the  electron 
beam  fluorescence  technique  indicated  the  possibility  of  rotational  non- 
equilibrium.In  these  studies  rotational  temperature  was  determined 
from  measured  rotational  fine  structure  of  the  emission  from  the  (0-0)  band 
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of  the  N2  first  negative  systeir..  They  all  determined  that  the  rotational 
temperature  deduced  from  rotational  line  intensities  was  higher  than  the 
calculated  or  inferred  translational  temperature  (inferred  from  auxilliary 
measurements  of  pitot  pressure  and  stagnation  conditions).  Furthermore, 
the  extent  of  apparent  nonequilibrium  between  the  rotational  and  trans- 
lational modes  (1)  increased  as  static  temperature  decreased  and  (2)  increased 
with  increasing  rotational  quantum  number,  i.e.,  fairly  close  agreement  be- 
tween rotational  and  translational  temperature  existed  down  to  20-30  °K  if 
only  the  first  four  rotational  line  intensities  were  used,  but  a higher  tem- 
perature resulted  from  the  first  eight  line  intensities,  and  an  even  higher 
one  if  additional  rotational  lines  were  incorporated.  Thus,  it  appeared  that 
rotational  temperature  increased  with  increasing  rotational  quantum  number. 

The  primary  emphasis  prevailing  at  that  time  was  that  this  anomaly 

o r 

was  the  result  of  either  (1)  inadequacies  of  the  Muntz  model  in  describing 

the  spectra  derived  from  electronic  excitation-deexcitation  processes'^  ’ 

or  (2)  the  presence  of  secondary  electrons  which  preferentially  excited  the 
33  34 

higher  J states,  ’ or  both.  While  evidence  was  presented  that  both  causes 

had  validity,  the  results  were  not  entirely  conclusive  that  some  rotational 

31 

nonequi’’ ibrium  did  not  exist,  except  for  Ashkenas  whose  tests  were  performed 

34 

in  an  essentially  static  gas.  Marrone  did  suggest  the  possibility  that  the 

apparent  rotational  nonequilibrium  might  be  real  for  his  free  jet  flow. 

Ashkenas  was  limited  to  liquid  nitrogen  temperatures  so  the  degree  of  non- 

^ r 0 1 

equilibrium  observed  was  smaller  1.1)  than  that  derived  in  studies 

trans 

performed  in  free  jets  or  nozzle  expansions  where  translational  temperatures 

down  to  4°K^^  were  inferred  - = 1.5).  It  is  possible  then  that  some 

trans 

fraction  of  the  apparent  discrepancy  between  the  rotational  and  translational 
temperatures  arose  from  the  electron  beam  excitation-deexcitation  model  and 


39 


the  rest,  particularly  at  the  lower  temperatures,  was  real. 

Vibrational  populations  were  computed  from  the  measured  gains 
using  a variety  of  rotational  temperature  models  which  had  the  general 
characteristic  that  T increased  gradually  as  J increased.  The  variation  of 
T that  appeared  to  produce  the  least  scatter  is  given  in  table  III.B.l. 

The  convention  used  here  is  that  J corresponds  to  the  rotational  quantum 
number  of  the  final  state,  e.g.  T(J=5)  is  the  assumed  temperature  corres- 
ponding to  the  relative  populations  of  rotational  states  J = 4 and  J = 5. 

It  was  necessary  to  assume  temperatures  T(l)  - T(4)  in  order  to  compute 
the  rotational  partition  function  (Equation  III.B.3).  The  contributions 
to  from  states  J ^ 11  are  less  than  .02%,  and  gains  were  measured  only 
for  5 £ J £ 10,  so  no  model  for  higher  level  temperatures  is  warranted. 

TABLE  III.B.l 

J 1-4  5 6 7 8 9 10 

T(J)°K  30.3  30.7  31.8  33.3  34.7  36.1  37.7 

The  results  using  this  model  are  shown  in  Figure  16.  The  scatter  has 
been  reduced  to  the  extent  that  all  points  fall  within  + 10;i  of  the  centroid 
of  values  computed  for  a particular  n,^  with  only  two  exceptions.  The  data 
fall  so  close  to  one  another  that  the  symbol  designation  of  Figure  15  becomes 
meaningless  in  Figure  16.  The  two  exceptions  are  that  values  of  n^  computed 
from  V=7  band  data  fall  a factor  of  two  (versus  a factor  of  four  in  Figure 
15)  below  those  computed  from  V = 6 band  data.  And  two  of  the  values  of 
Og  computed  from  V = 9 band  data  fall  distinctly  below  those  derived  from 
V = 8 band  results.  The  values  of  n^  computed  from  V = 6 and  9 band  data 
may  be  less  reliable  than  those  from  V = 7 and  8 bands  since  the  gains  are 
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smaller  and  thus  more  subject  to  statistical  error.  That  is,  solution  of 

Equation  III.B.l  for  a pair  of  gain  values  corresponds  to  seeking  a small 

difference  between  two  large  numbers  to  determine  n^  and  n,^  .j . Therefore 

small  variations  in  measured  gains  can  produce  large  changes  in  the  computed 

populations.  (Likewise  small  variations  in  assumed  rotational  temperature 

can  produce  large  variations  in  computed  populations.)  The  scatter  in 
28 

the  measured  gains  for  the  smallest  gains  is  sufficient  to  account  for  the 

apparent  disagreement  in  values  of  n^  and  ng  noted  above. 

The  high  degree  of  correlation  between  pooulations  computed  using 

a variety  of  rotational  levels,  we  believe,  justifies  the  assumption  of 

rotational  nonequilibrium.  There  is  no  independent  method  to  determine  if 

the  values  of  T in  Table  III.B.l  are  the  correct  ones  from  the  present  data. 

They  do  increase  gradually  from  T^^g^s  ’ would  expect,  and  they 

are  consistent  with  electron  beam  derived  results  for  similar  gas  tempera- 
30  33 

tures  ’ in  pure  nitrogen  flows.  Thus  we  conclude  that  when  the  rotational 
spacing  becomes  of  the  order  of  or  greater  than  the  translational  temperature 
in  supersonic  flows,  rotational  nonequilibrium  will  be  observed  even  when 
traditional  collision  frequency  and  rotational  equilibrium  concepts  indicates 
that  it  is  unlikely.  In  the  present  case,  for  example,  each  CO  molecule 
suffers  approximately  10  collisions/cm  transiting  the  test  section  which  is 
more  than  20  cm  long.  Several  measurements  in  shock  structure  studies  indi- 
cate that  two  tc  ten  collisions  are  required  for  rotational  equilibrium. 
However,  this  result  is  probably  misleading  in  the  present  case  since  in  a 
shock  wave  the  translational  temperature  exceeds  the  rotational  temperature 
in  contrast  to  the  rapid  expansion  ,,ase.  Moreover,  on  the  high  temperature 
side  of  the  shock, detai led  measurements  nave  not  been  performed  at  high 
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rotational  levels  where  the  rotational  spacing  2Je  might  be  greater  than 
the  translational  temperature  to  see  whether  they  are  indeed  equilibrating 
as  quickly.  Further  evidence  of  the  difference  between  shock  wave  and 
expansion  flow  rotational  equilibrium  has  been  recently  described  for  hydro- 
gen.  The  results  of  the  analysis  in  Reference  36  are  consistent  quali- 
tatively with  the  phenomena  suggested  here. 

It  is  interesting  to  note  that  the  rotational  nonequilibrium  ob- 
served in  highly  expanded  flows  presents  a serious  limitation  on  the 
achievable  small  signal  gain  or  output  power  that  can  be  derived,  'hat  is, 
while  we  certainly  observed  gains  that  were  large  compared  to  those  anti- 
cipated at  higher  temperatures  (less  expansion  from  identical  plenum  con- 
ditions), they  are  not  as  large  as  those  possible  if  rotational  equili- 
brium with  the  translatT  onal  mode  were  maintained. 


III.C.  Waveguide  Laser 


1.  Design  and  Operation 

A transverse  discharge  waveguide  laser  sketched  in  Figure  17  was 
constructed  using  a liquid  nitrogen  cooled  copper  block  as  one  electrode, 
a polished  aluminum  coated  glass  plate  as  the  other  electrode  and  edge 
polished  glass  slides  as  spacers.  The  waveguide  channel  was  2 mm  x 2mm 
square  and  10  cm  long.  Operation  was  achieved  at  a nominal  gas  tempera- 
ture of  100  °K  with  gas  mixtures  of  He: CO  in  the  range  of  4:1  to  10:1  at 
total  flow  rates  of  about  4 liters/min.  Upstream  and  downstream  pressures 
of  65  Torr  and  25  Torr,  respectively,  were  maintained.  This  corresponds 
to  a gas  exchange  rate  of  about  1.5  changes/  sec.  The  gas  discharge  was 
achieved  using  500  volt  pulses  which  had  a half  width  of  60  usee. 

The  laser  output  was  measured  to  peak  at  about  5-10  mW  with 
a Pulse  width  of  50  -120  usee.  The  electrical  efficiency  of  this  system 
based  on  an  assumed  peak  current  value  cf  4.6  amps  is  approximately 
.OOn.  It  was  not  possible  to  optimize  the  output  parameters  with  res- 
pect to  gas  mixture  or  discharge  conditions  because  of  the  deterioration 
of  the  strip  electrode. 

2.  Gas  Impurities 

During  the  course  of  the  gain  measurement  it  was  observed  that  a 
visible  layer  of  material  was  deposited  on  all  of  the  waveguide  walls 
if  the  discharge  was  operated  for  a half  hour  or  more.  For  a discharge 
period  of  about  30  minutes  the  deposits  appeared  as  a thin-film,  yellow 
and  red  in  color.  As  the  discharge  times  became  longer  the  deposited 
film  became  thicker,  more  uniform  and  brown  in  color.  One  of  these 
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deposits  is  pictures  in  Figure  18.  The  deposits  did  not  appear  in  a 
pure  helium  discharge  but  only  when  carbon  monoxide  was  present.  They 
also  did  not  appear  unless  the  discharge  was  switched  on;  i.e.  no 
deposits  were  left  from  the  flowing  gas  alone.  The  character  of  deposits 
was  also  not  affected  by  the  discharge  polarity. 

These  observations  suggested  that  there  was  some  impurity  present 
in  the  carbon  monoxide.  Since  the  same  bottle  of  CO  had  been  used  through- 
out the  research  period  it  was  conjectured  that  the  metal  walls  inside  of 
the  gas  bottle  could  have  been  'flaking'  and  progessively  adding  a metal 
(most  likely  iron)  impurity  to  the  gas.  (It  should  be  noted  that  a dry 
ice/methanol  cold  trap  was  introduced  in  the  gas  flow  line  to  facilitate 
the  removal  of  the  impurity.  By  visual  observation  the  character  of  the 
discolorations  was  unaffected  by  this  step.) 

The  possible  existence  of  an  impurity  was  investigated  by  examin- 
ing two  anode  plates  via  x-ray  spectrum  measurements  in  a scanning  elec- 
tron microscope  (SEM).  The  first  plate  (fl ) was  used  in  a discharge  for 
a total  of  about  12  hours.  The  initial  bottle  (A)  of  CO  was  used  to 
supply  this  discharge.  Thick  deposits,  brown/black  in  color  were  present 
on  the  plate  when  it  was  removed.  The  second  plate  (=2)  was  used  for 
about  2 hours  with  a new  bottle  (B)  of  CO  of  nominal  purity  99.99".  No 
colored  deposits  were  present,  however  a thin  whitish  film  and  bright 
spots  along  the  waveguide  path  were  observable.  For  an  operating  time 
of  two  hours  with  the  initial  bottle  of  CO,  colored  films  were  clearly 


present.  By  this  observation,  then,  the  purity  of  bottle  B was  much 
improved  over  that  of  bottle  A. 
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The  SEM  data  are  shown  in  Figures  19  through  24.  Referring  to 
Figure  18  there  were  three  regions  that  were  examined:  the  anode  plate 

off  of  the  deposited  stripe,  the  background  of  the  stripe  and  some  bright 
spots  which  appear  along  the  stripe.  For  both  plates  the  region  off  the 
stripe  showed  aluminum  only  (Figures  19  and  22).  In  the  case  of  plate  ?1 
the  stripe  background  is  clearly  composed  of  iron  (Figure  20)  while  the 
bright  soots  are  predominantly  copper  with  a substantial  amount  of  iron 
also  present  (Figure  21).  Both  regions  also  shown  the  presence  of  aluminum, 
as  expected.  The  data  for  plate  #1  show  that  the  signal  from  iron  is 
about  as  strong  as  that  for  aluminum  on  the  stripe  region.  (The  scanning 

o 

electrons  penetrate  several  hundred  A).  For  plate  ?2  the  strip  background 
contains  a slight  amount  of  iron,  approximately  .1%  of  the  amount  of 
aluminum  (Figure  23).  The  bright  spots  on  plate  #2  do  not  contain  any 
copper,  but  do  contain  an  increased  fraction  of  iron,  about  6%  (Figure 
24). 

From  these  data  we  infer  the  following: 

a)  Iron  is  present  as  an  impurity  in  both  bottles  of  carbon 

07 

monoxide  perhaps  in  the  form  FE(CO)0.  The  amount  of  the  impurity  present 
in  the  gas  probably  increases  with  time  because  of  'flaking'  of  the  walls 
of  the  gas  bottle.  Bottle  A was  in  use  for  one  and  one-half  years,  the 
impurity  became  noticeable  after  about  6 to  8 months.  In  addition,  the 
nominal  purity  of  bottle  A was  99.00%.  Thus,  the  reduced  fraction  of 
ircn  in  bottle  B versus  bottle  A is  most  likely  a consequence  of  its  higher 
purity  and  ' newness ' . 

b)  The  presence  or  absence  of  copper  in  the  bright  spots  probably 
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reflects  the  length  of  time  the  discharge  was  operating  and  the  severity 

of  arcing  in  the  discharge.  Cathode  sputtering  is  a common  phenomenon 

in  gas  discharges  wherein  the  impact  of  ions  on  the  catlicde  surface  causes 

some  cathode  material  to  be  removed.  The  high  current  densities  associated  ] 

with  arcing  enhance  the  sputtering  effect.  (Sputtering  at  the  anode  is 

not  likely  because  electrons  are  much  less  massive  than  the  positive  ions). 

Since  in  this  work  the  copper  block  serves  as  cathode,  sputtering  is  the 
likely  source  of  copper  on  plate  #1. 

The  absence  of  copper  on  plate  #2  reflects  the  different  con- 
ditions of  plates  #1  and  #2.  Plate  ^1  was  used  on  three  or  four  con- 
sectuve  occasions  for  a cumulative  discharge  time  of  about  12  hours.  During 
the  last  occasion,  arcing  was  especially  acute  as  a result  of  the  thick 
deposits  on  the  anode  and  cathode.  In  comparison,  plate  42  was  used  for  two 
hours  in  a very  homogeneous  discharge.  Arcing  was  present  but  greatly 
reduced. 

3.  Gain  Measurements 

Initial  measurements  of  the  gain  of  several  vibrational-rotational 
transitions  for  a gas  temperature  of  215  °K  exhibited  uniformly  small 
values  with  a peak  of  .064?i/cm  + .005  on  the  8 P(10)  line.  The  measured 
gains  were  generally  a factor  of  10  or  20  less  than  what  was  expected. 

These  magnitudes  did  not  vary  more  than  a factor  of  two  (i.e.,  .032% 
to  .C64?0  by  optimizing  gas  mixture  or  electric  discharge  properties. 

It  was,  however,  possible  to  verify  some  of  the  expected  qualita- 


tive features  of  the  small-signal  gain.  In  Figure  25  the  gain  is  plotted 
vs.  vibrational  level  for  three  sets  of  measurements.  The  magnitudes 
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are  normalized  with  the  peak  value  set  equal  to  one  in  each  case.  Two 
measurements  using  CO  gas  with  an  iron  impurity  (bottle  A)  both  shown 
peak  gains  in  the  range  V = 8 to  V = 10  and  a clear  reduction  at  V = 11 
and  V = 12.  The  higher  purity  CO  (bottle  B)  exhibited  a similar  dis- 
tribution of  gain,  but  with  the  peak  value  clearly  shifted  to  V=ll.  It 
should  be  recalled  that  previous  reports  of  gain  vs.  vibrational  level 
gave  V = 7,  8 and  V = 8,  9 as  the  optimum  values  of  V. 

The  distribution  of  gain  vs.  J-levels  is  not  so  clearly  defined. 
Measurements  were  made  on  the  8th  vibrational  level  for  J = 8,  9,  10,  12. 
Only  these  J-levels  were  used  because  lower  J-levels  have  transition 
frequencies  which  are  near  those  of  the  upper  J-levels  of  the  V = 7 
the  the  line  selection  mechanism  of  the  probe  laser  did  not  allow  clear 
resolution  between  these  transitions.  The  data  are  shown  as  total  per- 
centage gain  vs.  J,  in  Figure  26.  To  within  experimental  error  there  is 
no  apparent  gain  variation  over  J = 8,  9,  10  and  there  is  an  apparent 
decrease  at  J = 12. 


Increasing  the  O2  'f'low  rate  from  zero  to  about  20%  of  the  CO 

flow  rate  was  found  to  enhance  the  measured  gain  by  approximately  a factor 

of  two.  This  is  shown  in  Figure  27,  where  total  gain  (%)  is  plotted 

against  (liters/min).  Because  other  researchers  have  found  laser 
^2 

^ 38 

operation  to  be  sensitive  to  the  O2  fraction  it  was  expected  that  there 
would  be  an  optimum  value  of  with  gain  decreasing  on  either  side  of 

U2 

the  optimum.  This  is  shown  in  the  figure  by  a curve  extrapolated  to 


zero  gain  as  Qq  increases.  However,  the  data  points  at  Qq  /C^q  = 1.4 

and  2.7  indicate  that  the  expected  decrease  of  the  gain  as  Q-,  increases 

^2 
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does  not  occur.  Thus,  the  increase  of  gain  with  small  values  of  Q. 
agrees  with  observed  enhancement  of  CO  laser  oerformance  by  addition  of 
0,.  But  the  fact  that  the  gain  does  not  decrease  for  large  values  of 
Q,-,  is  contrary  to  expectations  and  not  well  understood. 

The  gas  impurity  discussed  in  the  previous  section  did  not  sig- 
nificantly affect  the  gain.  The  peak  total  gain  using  bottle  B was  .901^ 
whereas  using  bottle  A the  peak  was  .58%.  As  shown  in  Figure  25  the 
gain  peak  shifted  to  V = 11,  from  the  region  V = 8 to  V = 10  by  sub- 
stituting the  'clean'  CO. 

The  small-signal  gain  profile  of  the  waveguide  discharge  was 
found  to  be  spatially  non-uniform.  The  variations  occurred  only  along 
the  line  between  the  anode  and  cathode  (vertical  direction)  and  not  in 
the  horizontal  directions.  The  sequence  of  oscillographs  shown  in 
Figure  28  follows  the  movement  of  the  beam  center  across  the  waveguide 
from  anode  to  cathode.  As  the  beam  is  moved  in  the  direction  o'!"  the 
cathode  the  gain  increases  while  absorption  tends  to  zero;  when  the  beam 
moves  toward  the  anode  the  reverse  is  true.  This  was  observed  to  be  the 
case  on  several  transitions  and  with  the  discharge  polarity  reversed. 

Plots  of  measured  gain  and  absorption  versus  separation  between 
beam  center  and  waveguide  center  are  shown  in  Figures  29  and  30.  The 
measurements  of  Figure  29  were  made  for  an  average  pressure  of  44  Torr, 
those  of  Figure  30  for  82  Torr.  The  lOP  (10)  line  was  used  for  both 
measurements.  To  draw  conclusions  from  this  data  it  is  first  necessary 
to  recall  that  the  probe  beam  exhibits  a Gaussian  intensity  distribution 
as  shown  in  Figure  31.  Thus,  when  the  beam  is  aligned  for  maximum 
transmission,  the  peak  intensity  passes  through  the  waveguide  center 


while  at  the  walls  the  intensity  has  dropped  by  some  fraction. 

Consideration  of  Figures  29  and  30  shows  that  the  region  of  high 
gain  and  zero  absorption  appears  when  the  beam  center  is  just  past  the 
cathode  edge  of  the  waveguide  whereas  the  high  absorption  and  low  gain 
region  peaks  near  the  anode  wall  (1.0  mm  from  center)  and  tends  to  zero 
on  either  side  of  this  point.  This  indicates:  1)  that  the  region  of 
absorption  extends  over  a larger  portion  of  the  discharge  than  does  the 
gain  region;  ii)  that  the  gain  region  is  located  very  near  the  cathode. 
This  conclusion  is  drawn  because  the  gaussian  intensity  distribution 
drops  from  one  waveguide  edge  to  the  other  when  the  beam  center  is  posi- 
tioned outside  of  the  waveguide.  Consequently,  the  intensity  near  the 
cathode  wall  is  sufficient  to  exhibit  the  gain  region  but  decreases 
away  from  the  cathode  so  that  the  absorption  signal  is  very  small.  Fur- 
thermore, the  gain  goes  to  zero  if  the  beam  center  is  located  in  the 
waveguide  at  position  Z =-.75  mm.  Then  the  signal  intensity  is  smaller 
near  the  cathode  than  in  the  absorption  region. 

Two  further  points  should  be  noted.  First  of  all,  increasing  the 
pressure  appears  to  reduce  the  extent  of  the  gain  region.  Secondly, 
there  is  a region  of  absorption  and  gain  which  is  observed  after  the  beam 
center  passes  the  anode  edge  of  the  waveguide.  This  may  be  the  result 
of  grazing  reflections  on  the  anode. 

The  reason  for  the  non-uniformity  of  the  gain  is  not  certain.  One 

possible  cause  is  the  input  electrical  pulse  length.  The  gain  measurements 

were  performed  using  pulse  lengths  on  the  order  of  1 usec.  Computations 
0 

by  Boulnois  show  that  under  constant  pumping,  steady-state  vibrational 
populations  are  achieved  after  about  40  usec.  Furthermore,  it  should 
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be  recalled  that  laser  oscillation  was  achieved  using  100  usee  wide  pulses. 
It  is  possible  that  the  short  pulses  (1  psec)  do  not  deposit  enough  energy 
into  the  medium  for  inversion  to  be  obtained.  The  presence  of  gain  in 
the  cathode  fall  region  may  result  from  the  comparatively  high  electron 
energies  found  there.  That  is,  there  may  be  a spatial  distribution  of 
available  electrical  energy  such  that  gain  is  obtained  in  the  more  ener- 
getic cathode  fall  region  but  not  in  the  less  energetic  positive  column. 

On  the  other  hand  this  hypothesis  does  not  account  for  two  facts. 
First  of  all,  in  the  gain  region  the  absorption  is  zero.  The  expected 
mechanism  of  population  inversion  involves  an  initially  non-inverted 
(absorbing)  distribution  created  by  electron  energy  transfer  to  CO, 
followed  by  V-V  pumping  to  achieve  the  inverted  population.  Secondly, 

transverse  discharge  CO  lasers  have  been  reported  which  operate  using 

39  40 

voltage  pulse  widths  on  the  order  of  0.5  ysec.  ’ 
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I.  INTRODUCTION 

In  the  present  paper  we  shall  report  some  quantitative  information  on  CO-He 
VT  rates  at  high  vibrational  quantum  numbers  (i  = 25)  inferred  from  a collection 
of  quasi-equilibrium  CO  vibrational  population  measurements  recently  obtained  by 
Lightman  and  Fisher^  at  150°K. 

High  quantum  number  rates  are  difficult  to  measure  experimentally.  Know- 
ledge of  such  rates  are  important  in  kinetic  calculations  of  CO  laser  systems. 

In  the  absence  of  reliable  data,  they  are  usually  extrapolated  from  the  relatively 

abi-ndant  low  quantum  number  rates^’^’^’^’^’^’^’^’^^  using  theoretically  motivated 
11  12  13 

scaling  laws.  ’ ’ It  is  known  that  the  location  of  the  so-called  Plateau- 

14  15 

Boltzmann  knee  is  strongly  affected  by  the  high  quantum  number  rates  assumed.  ’ 

In  the  present  paper,  we  shall  take  advantage  of  this  sensitivity  to  infer  some 
quantitative  information  on  these  rates  from  experimentally  measured  quasi-equi- 
librium CO  population  distributions.  As  shall  be  shown  later,  this  indirect 
method  can  yield  a certain  ratio  of  VT  to  VV  rates  at  high  quantum  numbers.  How- 
ever, because  of  the  limited  amount  of  experimental  data  available,  we  are  able 
to  infer  information  only  near  i = 25  at  T = 150°K.  Nevertheless,  the  method  of 
data  reduction  is  interesting,  and  the  results  obtained,  although  limited,  appear 
to  be  the  only  really  high  quantum  number  (i  > 20)  experimental  CO-He  VT  data 
available  at  this  time. 

Recently,  Brechignac^^  using  a double  resonance  technique  made  measurements 
of  CO-CO  VV  rates  at  high  quantum  numbers.  One  of  his  major  conclusions  was  that 
the  near-resonant  VV  rates  extrapolated  to  high  quantum  numbers  by  the  Sharma- 
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12 

Brau  model  appeared  to  be  too  high.  For  CO-Ke  VT  rates,  the  9 i <.13  data 

3 

by  Hancock  and  Smith  appears  to  be  the  only  moderately  high  quantum  number 
data  available  in  the  literature  at  this  time.  A common  VT  rate  extrapolation 
scheme  suggested  by  Rich  et  al.^^  utilized  the  well-known  SSH^^  theory.  Recently, 
Verter  and  Rabitz  presented  detailed  theoretical  (non-linear)  calculations  of 
CO-He  rates,  and  their  results  at  high  quantum  number  and  low  temperatures 
(T  < 300°K)  disagreed  substantially  from  the  Rich  extrapolation.  In  particular, 
the  low  temperature  Verter-Rabi tz  VT  rates  increase  with  i much  more  rapidly  at 
low  quantum  numbers  and  much  more  gradually  at  high  quantum  numbers  than  that 
predicted  by  the  Rich  extrapolation.  The  present  work  provides  some  limited 
evidence  of  the  qualitative  features  of  the  Verter-Rabi tz  rates. 

II.  THE  EXPERIMENTAL  DATA 

The  Lightman  and  Fisher  data^  of  CO  vibrational  population  distributions 
mreasured  in  a CO-He  electric  discharge  is  reproduced  in  Figs.  1.  The  temper- 
ature of  the  system  is  145-1 55°K.  For  details  of  the  experimental  procedures 
the  readers  are  referred  to  the  original  paper. 

The  7-torr  and  4-torr  He  data  are  presented  separately  in  Figs,  la  and  lb, 

A 

respectf vely.  It  is  seen  that  the  "Plateau"  begins  at  approximately  i = 10, 
and  the  "Plateau-Boltzmann  knee"  occurs  near  or  beyond  i = 30,  the  upper  limit 
of  the  data.  The  7-torr  data  includes  various  combinations  of  CO  pressures  (.5, 
10,  1.5-torr)  and  discharge  power  inputs. 

Lightman  and  Fisher^  also  performed  detailed  computations  of  the  CO  popula- 
tion distributions  using  a computer  code.  The  CO-CO  VV  rates  assumed  were  based 
on  the  Sharma-Brau  scaling,  and  various  CO-He  VT  rates  were  used.  They  found 
that  neither  the  extrapolation  formula  suggested  by  Rich  et  al.  nor  the  theore- 
tical rates  by  Verter  and  Rabitz  could  generate  acceptable  agreements  with  the 
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measured  data;  the  former  produced  a Plateau-Bol tzmann  knee  at  too  low  values  of 
i while  the  latter  at  too  high  values  of  i.  When  the  Verter-Rabi tz  rates  were 
augmented  empirically  by  a multiplicative  factor  (i+1)^’^^,  excellent  agreements 
with  the  experimental  data  were  obtained.  (See  Fig.  lb)  Note  that  this  augmen- 
tation factor  is  approximately  8 when  i = 25. 
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f 
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is  obtained  from  the  dipole-dipole  transition  matrix  elements  computed  by 

Young  and  Eachus.^^  It  can  be  shown  that  so  computed  may  adequately  be 

VV 

represented  by  taking  C =1.5. 

VT 

The  VT  adiabaticity  factor  is  generally  expected  to  be  an  exponen- 
tially increasing  function  of  i.  Normalizing  so  that  A^^  = 1 , we  can  de- 
termine u)'^'^(150°K)  = 1.5  X 10  cm Vsec-molecule  from  the  available  experi- 
2 

mental  data.  Although  the  Verter-Rabitz  rates  are  computed  from  a fully 

1 9 

non-linear  theory,  their  rates  can  be  correlated  by  Eq.  (3.1).  At  150°K, 
the  correlation  yields 

A^'*'(150°K)  = a^-  exp(0.09i),  (Verter-Rabitz^^),  (3.5) 

where  a^(150°K)  is  a rapidly  increasing  function  of  i in  the  first  few  levels 
and  tends  to  a constant  value  (a^.  (150°K)  = 4.5)  for  i 2.  10.  According  to  the 
Rich  extrapolation,  A^^  (150°K)  is  given  by 

Ay^(150=K)  = exp(0.24i),  (Rich  et  al.^^  (3.6) 


These  two  rates  are  seen  to  disagree  strongly  with  each  other  at  high  quantum 
numbers . 

Since  the  forward  and  the  reverse  rates  are  related  by  detailed  balance. 


A^'|(T)  is  often  rewritten  as 

‘ J 


A^j(T)  = P^^(T)  exp 


(AEj-AE.)/2kT 


(3.7) 


,VV 


th 


where  P]j(T)  is  a symmetric  matrix,  E^  is  the  vibrational  energy  of  the  i 
level  (E.  = i(l-a(i-l))k0.^.^,  = 3084°K),  and  AE.  = E.^.,  - E.  = (1 -2ai  )k9^.^. 

Normalizing  P^j(T)  so  that  P^^(T)  = 1 , we  can  determine  1 50°K)  = 1.5  x 10 
cm^/sec-molecule  from  the  available  low  quantum  number  data.  We  shall  define 
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the  second  moment  of  P^'|(T)  by 

* V 


bi(T)  = j ^ 

j=0 

Theoretically,  P^j  should  depend  primarily  on  |i-j|  and  should  decrease 


(3.3) 


Thus,  b^(T)  should  be  insensitive  to  i in  the  plateau 


rapidly  wi th  | i-j 

region  and  shall  be  written  hereafter  simply  as  b(T).  Using  the  experimental 

and  theoretical  data  on  . (1 50°K)  collected  by  Brechignac,  we  obtain 

^ > J 


b(T)  = 6.6 


(l50  ) 


3/2 


(3.9) 


IV.  THE  THEORETICAL  PLATEAU  DISTRIBUTION 


21 


Using  a differential  approximation  first  introduced  by  Brau  and  later 

i 


14  15 

refined  by  Gordiets  et  al.  and  Lam,  the  Plateau  distribution  n.  of  CO  can 


be  explicitly  expressed  in  terms  of  the  assumed  kinetic  rates: 

VT  i 


T oVV 

^ = TV  IvV 

^i 


To)' 


He  VV 
^bgcu 


(4.1) 


where  r is  a reference  quantum  number,  n^^  is  the  number  density  of  He  which 


is  assumed  to  dominate  the  VT  processes,  and  g = AE^.  - = 3 

14  15  V*^ 

Eq.  (4.1)  differs  slightly  from  earlier  works  ’ by  recognizing  that  ' 

and  may  be  different  as  mentioned  earlier. 

We  shall  represent  at  high  quantum  numbers  by 


7°K. 


A^  (150°K)  = B^  exp  S(i-r) 


(4.2) 


where  B^  and  3 are  constants  to  be  determined.  Since  1^/1^  is  only  weakly 
dependent  on  j , we  shall  approximate  it  by  • Eq.  (4.1)  can  now  be  re- 
written as 


n . = 
1 


1 

,VV 


nVV  P 

"r*r  - "He 


jSU-r)., 


1 -e 


-6 


(4.3) 
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where 


0 T VT  .VT 
T u)  Jl 

Zbg  oj 

Eq.  (4.3)  contains  only  three  kinetic  rate  related  constants,  , C,  and  3. 


(4.4) 


V.  DATA  REDUCTION 

The  data  reduction  procedure  is  as  follows.  From  a given  measured  distri- 
bution, the  value  of  is  first  determined  to  fit  the  data  at  i = r,  and  then 
the  values  of  C,  , and  3 are  determined  by  trial  and  error  to  obtain  accept- 
able agreement  over  the  range  10  <_  i ^ 30.  Choosing  r = 25,  it  is  found  that 

-4 

so  long  as  C is  within  15%  of  1 .2  X 10  , acceptable  fits  can  be  generated  with 

= -1  to  +1.5  and  3 = 0.001  to  0.11.  In  other  words,  within  the  ranges  indi- 

VV 

cated  Eq.  (4.3)  is  not  sensitive  to  the  values  of  C and  3,  but  is  sensitive 

mainly  to  the  value  of  C.  The  CO  population  distributions  calculated  by  Eq.  (4.3) 
-4  VV 

using  C = 1.2  x lo  , C =1.5,  and  3 = 0.09  are  shown  as  dotted  lines  in  Figs.  1. 
The  solid  line  in  Fig.  lb  is  the  exact  numerical  solution  obtain  by  Lightman  and 
Fisher  using  the  Verter-Rabitz  rates  augmented  by  (i+1)^’^^  as  mentioned  earlier. 

V.  RESULTS  AND  DISCUSSION 

From  Eqs.  (3.1),  (4.2),  and  (4.4)  the  CO-He  VT  rates  at  high  quantum  num- 
bers can  be  written  as 

VV  „VV 

VT  ^ VT 

k^^(if)  = wj C exp  B(i-r)  . (5.1) 

T Ji;'  ' 

It  is  seen  that  for  i r = 25,  the  value  of  k'^^(i-l-)  is  relatively  insensitive 

VV 

to  the  uncertain  values  of  C and  3 in  the  indicated  range.  Using  g = 37°K 
and  the  150®K  data  for  b,  C quoted  earlier,  we  have 
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^VV  • 

k''"'’(U;150°K)cmVsec-niolecule  = 6.8  x exp  6(i-25) 


vT  n 

‘'25 


(5.2) 


vv 

Usi.'ig  C = 1.5,  we  have 


k'^^(i  + ;150°K)  = 1.3  x -Jr  S(i-25) 

ft  * * 

^25 


(5.3) 


which  is  the  150°K  CO-He  VT  rates  inferred  from  the  Lightman  and  Fisher  data. 

The  Verter-Rabitz  and  the  Rich  extrapolation  rates  can  be  written  in 
similar  forms: 


Verter-Rabitz 

k''^(i;;150°K) 
Rich  Extrapolation 


k''^(if  ;150°K) 


= 1.8  X 


exp 


0.24(i-25) 


(5.4) 


(5.5) 


Thus  it  is  seen  that  near  i = 25,  the  Verter-Rabitz  rates  are  lower  in  magni- 
tude than  the  present  experimentally  inferred  rates  by  a factor  of  about  7, 
but  their  slow  variation  with  i (B(Verter-Rabitz)  = 0.09)  is  consistent  with 
the  Lightman  and  Fisher  data.  On  the  other  hand,  the  Rich  extrapolation  rates 
intersect  the  experimentally  inferred  rates  near  i = 25,  but  deviates  sub- 
stantially from  it  for  other  values  of  i because  of  the  relatively  large  value 

of  0 = 0.24  used.  In  the  Rich  extrapolation,  0 is  estimated  by  the  Keck  and 
1 3 

carrier  curve-fit  representation  of  the  VT  adiabaticity  factor: 


6 = I cte^-b  (5.6) 

where  = 3084°K  is  the  characteristic  vibrational  temperature  of  CO  and 
is  a constant  proportional  to  the  interaction  length  L of  the  short-range 
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epulsive  potential.  The  value  3(150°K)  is  obtained^^  by  using  = 0.114 

which  corresponds  to  L = O.zA.  It  is  of  interest  to  observe  that  if  a (nonnal- 

19 

ized)  detailed-balance  factor  is  empirically  included,  we  would  have 

, (5 

yielding  S(150°K)  = 0.114  which  is  consistent  with  the  present  results. 

Since  experimental  CO-He  1-0  VT  rates  are  wel 1 -documented,  the  small  value 

of  p and  the  relatively  large  magnitudes  of  rates  at  high  quantum  numbers 

VT 

inferred  from  the  Fisher-Lightman  data  suggest  that  must  increase  rapidly 

in  the  first  few  levels  (i  <10).  This  observation  is  qualitatively  consistent 

1 9 

with  the  Verter-Rabitz  results  as  correlated  by  Eq.  (3.5).  It  was  shown  that 
a^.  (T)  is  essentially  unity  for  T > 300°K,  but  exhibits  a rapid  i dependence  at 
the  first  few  levels  for  300°K  > T ^ 100°K.  This  interesting  and  unexpected 
behavior  of  the  Verter-Rabitz  results  has  so  f • r not  been  satisfactorily  inter- 
preted theoretically.  The  present  inferred  rates  provide  the  first  experimental 
indications  in  support  of  this  qualitative  result. 

It  is  important  to  note  that  the  value  of  C inferred  from  the  measured  data 
is  a ratio  of  VT  and  VV  rates.  Implicit  in  Eq.  (5.3)  is  the  assumption  that  the 
high  quantum  number  nearly  resonant  VV  rates  are  correctly  scaled  by  the  Young 
and  Eachus  transition  matrix  elements  (since  = 1.5  was  used).  Since  recent 

measurements  by  Brechignac^^  indicates  that  such  scaling  appears  to  overestimate 
the  high  quantum  number  VV  rates,  our  inferred  rates  will  be  similarly  affected. 

The  data  was  originally  reduced  by  choosing  r = 10  as  the  reference  quan- 

VV 

turn  number.  Although  the  combinations  of  C,  C , and  3 fitting  the  data  appear 
disparate,  the  inferred  CO-He  VT  rates  from  each  set  of  parameters  intersect 
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each  other  near  i = 25,  indicating  that  only  the  local  VT  rates  are  important 
for  the  given  collection  of  data.  The  value  of  r = 25  was  then  chosen  so  that 
the  curve-fitting  became  sensitive  mainly  to  the  value  of  C alone.  Since  the 
present  inferred  rates  are  obtained  from  a limited  collection  of  data,  addi- 
tional experiments  with  a much  wider  range  of  external  excitation  strength 
and  He  pressures  are  clearly  necessary  to  collaborate  and  refine  the  present 
findings. 


F 
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FIGURE  CAPTIONS 

Fig.  la  CO  vibrational  population  distributions  measured  by  Lightman 

and  Fisher^  with  7-torr  He.  Dotted  lines  are  Plateau  distri- 
butions computed  from  Eq.  (4.3)  using  C = 1.2  x 10"^,  = 1.5, 

j and  S = 0.09. 

Fig.  lb  CO  vibrational  population  distributions  measured  by  Lightman 

and  Fisher^  with  4-torr  He.  Solid  line  is  exact  solution 
computed  by  using  tne  CO-He  VT  rates  by  Verter  and  Rabitz 
augmented  by  (i+l)^’^^.  Dotted  line  is  Plateau  distribution 
computed  from  Eq.  (4.3)  using  C = 1.2  y lu~^,  = 1.5,  and 

6 = 0.09. 
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ABSTRACT 

An  analytical  theory  of  vibrational  relaxation  is  presented  showing 
explicitly  the  effects  of  the  strength  of  the  pumping,  the  anharmonicity  of 
the  energy  level  spacings,  and  the  VV  and  VT  kinetic  rate  models.  The  quasi- 
equilibrium population  distribution,  the  loci  of  the  Treanor-Plateau  boundary 
and  the  Plateau-Boltzmann  knee,  the  nominal  magnitude  of  the  plateau,  the 
ground-level  vibrational  temperature,  and  the  rates  of  vibrational  quanta  and 
energy  dissipations  are  obtained.  A method  for  deducing  quantitative  informa- 
tion on  kinetic  rates  at  high  quantum  numbers  from  measured  population  distri- 
butions is  suggested. 


li 


A 


I 


102 


AN  ANALYTICAL  THEORY  OF  VIBRATIONAL  RELAXATION 
FOR  ANHARMONIC  MOLECULES  UNDER  STRONGLY-PUMPED  CONDITIONS 

S.  H.  Lam"'" 

I.  INTRODUCTION 

When  a diatomic  gas  at  temperature  T is  subjected  to  external  vibrational 
excitation  (such  as  an  electric  discharge),  its  vibrational  population  distri- 
bution N^  will  deviate  from  the  equilibrium  Botlzmann  distribution  given  by 

Ni^(T)  E N^  exp(-Ei/kT)  (1.1) 

where  E-  is  tne  energy  of  the  i^*^  vibrational  level.  The  non-equilibrium  N.^ 
is  governed  by  a system  of  master  kinetic  equations  and  depends  in  general  on 
the  strength  of  the  external  excitation,  the  anharmonici ty  of  the  energy  levels, 
and  the  relevant  kinetic  rates.  The  vibrational  relaxation  process  is  dominated 
by  the  fact  that,  in  some  average  sense,  VV  collision  rates  are  much  more  rapid 
than  VT  collision  rates.  Vibrational  quanta  externally  pumped  into  the  gas  can 
be  redistributed  by  numerous  quanta-conserving  VV  collisions  before  being  deac- 
tivated by  the  much  slower  VT  collisions.  We  shall  be  interested  in  the  quasi- 
equilibrium distributions  of  N.j  under  strongly  pumped  conditions. 


By  assuming  that  VT  rates  are  uniformly  negligible  in  comparison  to  VV 
rates,  Osipov^  showed  that  for  harmonic  oscillators  N.  would  remain  Boltzmann 
and  could  be  characterized  by  a single  vibrational  temperature  T^: 
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N^.  = N^.^(T^),  hannonic  oscillators.  (1-2) 

The  more  general  resu.t,  valid  for  arbitrary  energy  level  spacings,  ivas  ob- 

2 

tained  by  Treanor,  Rich,  and  Rehm  : 

N.  = N 5 N.^T)  exp(iy)  (1.3) 

where  the  constant  y replaced  as  the  single  parameter.  In  either  case,  a 

relatively  simple  relaxation  equation  could  be  derived  for  or  y. 

' These  simple  and  elegant  theories  serve  adequately  only  when  the  strength 

i of  the  external  excitation  is  weak.  Motivated  by  its  applications  to  molecular 

^ infrared  lasers,  the  complete  master  kinetic  equations  .have  rec-'ntly  been  ex- 

3 4 5 6 

tensively  studied  numerically.  ’ ’ ’ The  computed  under  strongly  pumped 

conditions  exhibits  the  following  general  features,  (a)  In  the  lower  levels, 

0 ^ i < i*,  indeed  can  be  approximated  by  ri^"''^^{T ;y)  for  some  y > 0.  (b)  In 

★ ★★ 

the  intermediate  levels,  i < i < i , is  relatively  flat  and  decreases 
slowly  with  i.  (c)  For  i = i , drops  sharply  and  tends  toward  (T). 

These  three  regions  shall  be  called  the  Treanor,  the  Plateau,  and  the  Boltzmann 
regions,  respectively.  These  qualitative  features  appear  to  be  relatively  in- 
sensitive to  the  kinetic  rate  models  used. 

3 

The  physical  interpretations  of  the  above  observations  are  clear. 

Normally,  vibrational  quanta  are  pumped  into  the  gas  mostly  at  low  levels. 

When  the  pumping  is  weak  and  the  input  quanta  can  be  deactivated  by  local  VT 
collisions,  the  Treanor  distribution  results.  When  the  pumping  is  sufficiently 
^ strong,  local  VT  deactivations  become  inadequate.  The  rapid,  quanta-conserving 

VV  collisions  aided  by  the  small  anharmonicity  of  the  energy  levels  then  prefer- 
entially transport  the  excess  quanta  up  the  vibrational  ladder  to  be  deactivated 
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by  the  much  faster  VT  colTisions  there.  The  Plateau  region  which  appears  only 

under  strongly  pumped  conditions  serves  primarily  as  a conveyer  belt  between 

the  source  and  the  sink  of  vibrational  quanta. 

In  the  present  paper,  we  present  a simple  theory  of  vibrational  relaxation 

for  anharmonic  oscillators  under  strongly  pumped  conditions.  In  contrast  to 

numerical  studies  where  kinetic  rates  and  other  input  data  must  be  explicitly 

specified,  our  analytical  results  are  expressed  in  terms  of  the  input  data, 

showing  clearly  the  effects  of  kinetic  rate  models.  The  analytical  techniques 

7 8 9 

used  follow  closely  the  works  of  Brau  and  Gordiets  et  al . ‘ The  quasi-equi- 
librium distribution  N^.,  the  loci  of  the  Treanor-Plateau  boundary  and  the 

* ★★ 

Plateau-Boltzmann  knee  N (i)  and  N (i),  the  ground-level  vibration  temperature 
T^^,  the  nominal  magnitude  of  the  plateau  N’^,  and  the  rates  of  vibrational  quanta 
and  energy  dissipation  D'^"^  and  are  obtained.  It  is  shown  that  only  a few 
parameters  suitably  extracted  from  the  large  collection  of  VV  and  VT  rates  are 
important.  A method  for  quantitative  determination  of  kinetic  rate  information 
at  large  quantum  numbers  from  experimental  N^.  measurements  is  suggested. 

II.  KINETIC  RATE  MODELS 

We  shall  assume  that  can  be  adequately  represented  by 

E.  = i l-a(i-l)  ke  ,,  (2. 

where  a is  the  anharmonicity  and  is  the  characteristic  vibrational  temper- 
ature of  the  diatomic  molecule  in  question.  We  shall  use  A to  denote  the 
forward  difference  operator.  The  energy  level  spacing  is  thus  given  by 

AE,  5 E,„  - E,  = 


(2.2) 
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Applying  A to  Eg.  (2.2),  we  have 
A^E.  E A(AE.)  = 

For  CO,  a = 5.98  x 10“^  and  0^^^  = 3034°K.  We  shall  see  that,  inspite  of  its 
smallness,  the  anharmonicity  a plays  an  important  role  in  the  later  develop- 


ments . 

The  major  kinetic  rates  of  interest  are  the  so-called  VV  and  VT  single- 
quantum-jump rates.  Let  k'^'^(it,j-l-)  denote  the  vibration-vibration  (VV)  rate 
of  a i -»■  i+1,  j+1  j collision  between  two  molecules  of  interest.  Let 
k^^(i+l  J i)  denote  the  vibration-translation  (VT)  rate  of  a i+1  t i collision 

a 

between  a molecule  of  interest  and  a partner  atom  or  molecule  identified  by 

the  subscript  a.  Guided  by  the  v;ell -known  SSH  theory, ^ we  shall  formally 
, VV  VT 

express  k and  k as  follows: 

d 

k''^if,j  + ) = J^T)^/j  exp  [(AEj-AE.)/2kT]  (2 

k^^(i+l  J i)  = o;^'^(T)^.  A^'[^  exp  [(±AE . -AE . )/2kT]  , (2 

where  Jt-,  A^^,  and  A^"*^  are  dimensionless  and  are  normalized  as  follows: 

1 1 j 1 ,a 

t = 1 a'^^  = 1 a'^^  =1  (2 

^0  '*  ^ii  ’’  %,a  ‘ • 

In  addition,  detai  led-balance  requires  = A^y. . The  dimension  of  k'^'^ 

1 J J 1 

and  k^^(cmVmolecule-sec)  is  carried  by  and  representina  k'^^(lt,  1 + ), 
k^"*^(l  0)  respectively. 

a 

Theoretically,  is  proportional  to  the  square  of  the  oscillator  matrix 
element  and  and  A^^,  are  the  adiabaticity  factors.  Physically,  gives 

IJ  1 id  1 

the  quantum  number  scaling  of  the  resonant-exchange  VV  rates: 


kVV(i^,  i;)  = JV 
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For  harmonic  oscillators,  can  be  shown  to  be  simply  i + 1.  We  shall  assume 
here  simply  that  is  known  and  that  is  small. 

The  VV  adiabaticity  factor  depends  in  principle  on  both  i and  j. 

* *J 

1 3 

Guided  by  the  theories  of  Landau  and  Teller  and  SSH,  we  shall  assume  that 


a"  . A^h-T) 

14 

where  C is  the  so-called  adiabasy  parameter  defined  by 


5 = oe,,b  |i-j|V 


e»''/T 


(2.8) 


(2.9) 


and  is  a constant  which  depends  on  the  reduced  mass  and  the  interaction 
potential  of  the  collision.  For  CO,  is  0(10  ^)°K  \ In  addition, 
we  shall  assume  that  a'^'^(5;T)  decreases  exponentially  with  increasing  C to  re- 
flect the  fact  that  nearly-resonant  collisions  are  much  more  rapid  than  non- 

VV 

resonant  collisions.  I'Jote  that  in  general  A^j  niay  also  depend  on  i + j but  liiis 
dependence  has  been  neglected. 

The  VT  adiabaticity  factor  is  also  expected  to  depend  on  an  analogous 

adiabasy  parameter.  We  shall  merely  require  that  Ay"*”,  increases  exponentially 

I ,a 

wi th  i . 

III.  THE  MASTER  KINETIC  EQUATIOHS 

We  consider  a diatomic  gas  with  number  density  n^  in  the  presence  of  an 
atomic  diluent  with  number  density  n^.  The  master  kinetic  equations  can 

a 

readily  be  written  in  the  following  form: 


dt 


L 


k=i+l 


pVV 


0,1,2,.. 


(3.1) 
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where  N.  is  the  normalized  vibrational  population  distribution  of  the  molecules 

CO 


Z N.  = 1,  (3.2) 

i=0  ^ 


VV  VT 

and  the  dissociation  level  is  represented  by  i = The  terms  , 
F^  represent  the  net  fluxes  of  m.olecules  (per  molecule  per  second)  going 
from  any  level  at  or  below  level  i to  any  level  at  or  above  i + 1 by  VV, 
VT  and  any  other  processes,  respectively.  The  VV  and  VT  fluxes  can  be 
written  as  follows; 


= "a 


k"(i  - - kf(HI  ^ 


1 


(3.3) 

(3.4) 


where  v;e  have  assumed,  for  the  sake  of  simplicity,  that  F^  is  dominated  by 
collisions  with  the  diluent  atoms.  In  general,  Fy"*"  should  be  the  sum  of  all 
VT  contributions  including  self-collisions. 

Using  the  VV  and  VT  rate  formulas  discussed  earlier,  we  can  rewrite 
Eqs.  (3.3)  and  (3.4)  as  follows: 


where  S..  is  a symmetric  matrix 

' J 


VV, 


Sij  = S.i.Si  A-(C;T)s., 


and  s^  and  f^  are  given  by 


s.  = ^.N.  exp  (-4E./2kT), 


fi  = s.  exp  (-AE./2kT) 


(3.5) 

(3.6) 


(3.7) 


(3.8) 


(3.9) 
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The  auxiliary  dependent  variable  is  defined  by 

_Nhi  N^T) 

' = "i  nB^,(T) 


= exp  + In  N.j. 

R TRR 

Note  that  = 1 and  M.  = exp  y correspond  to  N^  = N“(T)  and  (T;y). 

respectively. 


The  master  kinetic  equations  can  now  be  written  in  the  following  compact 


form: 


where 


I 'iJ^^-’'V  = ^''(M.-l)f.  -Q. 
j=0 

n VT 
VT  ^ ^a  ^a 

e - vv  ’ 

CO 

m 

Oi-T-TvtMt  t 

"m*"  k=i+l 


The  term  represents  mainly  external  excitation  (such  as  electric  discharge, 
radiative'  transitions,  etc.)  but  may  include  mul ti-quanta-jump  and  other  pro- 
cesses. The  dimensionless  Q.  combines  these  external  excitations  and  the  un- 
steady term  and  is  to  be  considerd  as  a forcing  term  in  Eq.  (3.11).  The  conven- 
tional master  kinetic  equations  can  be  recovered  from  Eq.  (3.11)  by  applying 
the  forward  difference  operator  A. 

Summing  Eq.  (3.11)  over  all  i,  we  obtain  the  so-called  quanta  conservation 
equation 


dt  Z.  M 


(3.1^) 
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where  a is  the  specific  quanta  (average  number  of  quanta  per  molecule)  defined 
by 

^ 00  00  oo 

” = I '"i ' Z Z "k  • ' 

i=0  i=0  k=i+l 

and  is  the  rate  of  quanta  deactivation  by  VT  collisions 


z vr  z • 


The  contribution  from  quanta-conserving  VV  collisions  vanishes  identically. 

Multiplying  Eq.  (3.11)  by  AE^  and  then  summing  over  all  i , we  obtain  the 
vibrational  energy  equation: 

00 

dT  ®vib  X ’"i'^^i  " *^vib  ’ 
i=0 

where  e^^j^  is  the  specific  vibrational  energy: 

CO  00  oo 

^(b  = Z ^i”i  ’ Z Z • 

i=0  i=0  k=i+l 

and  D ..  is  the  rate  of  dissipation  of  vibrational  energy  by  both  VV  and  VT 

VI  b 

collisions: 

i=0  i=0 


IV.  PRELIMItlARY  DISCUSSION 


We  may  consider  n and  n to  be  representative  frequencies  of  VV  and 
•'  m a a 

VT  collisions,  respectively.  The  dimensionless  parameter  can  thus  be 
interpreted  as  the  number  of  VT  collisions  per  VV  collision  and  is  in  general  a 


small  number.  We  shall  define 


'total 


^Z»i  = 


— y F=  - 

VV  1 dt 


n u)  ' . . 
m 1=0 


110 


which  may  be  loosely  interpreted  as  being  the  number  of  quanta  externally  pumped 
into  a molecule  of  interest  per  VV  collision.  Under  most  practical  situations 
of  interest,  Q^Q^-gi  is  also  a small  number.  We  shall  confine  our  attention  to 
strongly-pumped  conditions  defined  by 


1 >>  <^total  ^ 


VT 


(4.2) 


In  addition,  since  the  bulk  of  the  vibrational  population  resides  in  the  first 
few  levels,  we  shall  assume  that  occurs  mainly  at  low  levels.  Consequently, 
^total  approximated  by  summing  over  only  the  first  few  levels. 

If  we  formally  invoke  Eq.  (4.2)  and  completely  neglected  the  right-hand  side 
of  Eq.  (3.11),  we  would  obtain 


j=0 


(4.3) 


from  which  the  Treanor  solution  = constant  can  immediately  be  recovered  by 

inspection.  While  this  simple  and  elegant  solution  can  be  justified  for  weakly- 
VT 

pumped  (1  » e » Q^Q^g-|)  conditions,  it  is  known  to  be  incorrect  for  strongly- 
pumped  (1  >>  >>  conditions.  In  fact,  Brau^  has  shown  that  this  Treanor 

solution  is  not  unique  but  is  merely  a member  of  a family  of  solutions  for 
Eq.  (4.3). 

The  major  physical  simplification  underlying  the  present  development  is  the 
dominance  of  the  rapid  nearly-resonant  VV  collisions  over  the  slov.'er  non-resonant 
VV  collisions.  The  VV  adiabaticity  factor  a'^'^(5;T)  is  generally  expected  to  be 
some  exponentially  decreasing  function  of  5.  Consequently,  the  elements  of  the 


non-linear  symmetric  matrix  operator  S..  as  defined  by  Eq.  (3.7)  are  expected  to 

^ 0 

be  largest  in  the  vicinity  of  the  diagonal.  Under  the  assumption  that  a'’^'’^(;;T) 


decreases  sufficiently  rapidly  with  5 (see  Eqs.  (7.7)  later),  the  matrix  operator 


Ill 


S-.  can  be  converted  into  a differential  operator.  The  original  discrete  sys- 

* J 

tem  of  master  kinetic  equations,  Eq.  (3.11),  thus  becomes  a single  ordinary 
differential  equation  which  can  be  studied  in  details  by  standard  techniques.^ 

V.  THE  DIFFERENTIAL  APPROXIMATION 


,VV 


For  A (5;T)  decreasing  sufficiently  rapidly  with  5,  the  main  contributions 


-VV 


to  f!'  as  given  by  Eq.  (3.5)  will  come  from  nearly-resonant  (j  = i)  collisions. 


Accordingly,  we  formally  expand  M.  and  s-  about  j = i as  follows 

J W 


M,. 


1 


M.  + (j-i)AM.  + 2 (j-i)^A^M.  + ...  , 


Sj  = s,  ♦ (i-1)ASi  + ...  , 

Substituting  the  above  in  Eqs.  (3.5)  and  (3.7),  we  obtain  formally; 


(5.1) 

(5.2) 


rVV  _ „ VV 


ja,s;raH,  + b,A(s?AH.)  t ...  j 


(5.3) 


where 


ai(T)  H ^ (j-i)A'^^5;T) 

b^d)  H Y 1 (j-i)^A''^5;T)  . 
j=0 

For  values  of  i sufficiently  large  such  that  the  lower  limits  of  the  summation 
are  unimportant. 


i » ?(T)  E (T/e^^f^/ad 


vib  ’ 


(5.4) 


a^(T)  and  b^(T)  become  independent  of  i.  We  obtain 
a.(T)  = ^ nA'^^n/i;T)  = 0 , 


p = .CO  CO 


b(T)  = b.(T)  = ^ n^A''dn/i;T)  . 
n=0 


(5.5) 
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Thus  the  leading  approximation  to  can  be  written  as  follows: 

= bn^cu'^'^A(s^AM^ ),  i » i . (5.( 

/s 

Note  that  non-resonant  VV  contributions  (|j-i|  » i),  even  those  with  partners 
from  the  highly  populated  lower  levels,  have  been  neglected.  We  shall  re- 
examine this  approximation  later  in  Section  VII. 

Under  strongly  pumped  conditions,  the  number  of  levels  with  significant 
populations  can  be  relatively  large.  We  shall  take  advantage  of  this  observation 
and  represent  the  discrete  variables  , E^,  etc.  by  continuous  functions 

N(i).  M(i)»  E(i)»  etc.  We  further  assume  that 

A^E 

1 » |A  £n  N.|,  1 » 6 E - > 0 , (5.; 

so  that  A can  be  approximated  by  the  differential  operator  d/di.  For  example, 
as  defined  by  Eq.  (3.10)  can  be  represented  by 


M.  = M(i)  = exp  ^ + In  N 


di  Ut 


and  AM^  by 


Using  this  differential  approximation,  Eq.  (5.6)  can  be  written  as 


Note  that  the  VV  adiabaticity  factor  A^^(5;T)  is  represented  only  by  its 
second  moment  b(T)  which  can  be  written  alternatively  as: 


b(T)  = f y*  C=A^'^(r;T)dC  . 
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The  detailed  shape  of  is  seen  to  be  unimportant  so  long  as  it  decreases 

sufficiently  rapidly  with 
VT 

The  term  F^.  can  similarly  be  approximated  by  a continuous  function.  For 

strongly  pumped  conditions,  the  magnitude  of  is  generally  very  large  except 

★ ★ 

in  the  Boltzmann  region  where  it  is  approximately  1.0.  We  shall  denote  i to 
be  the  Plateau-Bol tzmann  boundary  defined  tentatively  by  M(i<i**)  >>  1 and 
M(i>i**)  =1-  tq.  (3.6)  can  thus  be  written  as 


" "a^^r  o < i < i** 

I a a a 


^0  . 


. ^ . ★★ 

■>  X”  T 

a. 


(5.12) 


Using  Eqs.  (5.10)  and  (5.12),  the  originally  discrete  master  kinetic  equations 
can  be  represented  by  the  following  differential  equation:^ 


r .2 

dTT 


= - Q,  i(T)«i<i**, 

a 


= 0 , 


(5.1 


where  Q(i)  is  the  continuous  representation  of  Q^. . 


VI.  THE  MOHiriAL  PLATEAU  DISTRIBUTION 

If  the  right-hand  side  of  Eq.  (5.13)  is  formally  neglected,  we  obtain 
upon  integration^ 


Inli  - 6 


= - C 


(6.1) 


where  C is  an  integration  constant.  In  spite  of  the  approximations  invoked, 
we  can  readily  verify  that  N(i;C=0)  is  precisely  the  Treanor  solution: 

N(i;C=0)  • =xp}i6(i- 


(6.2) 


(>-) 
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where  the  integration  constant  i^^^  is  the  location  of  the  so-called  Treanor 
minitnum.  It  is  clear  that  for  C / 0 a family  of  ncn-Treanor  solutions  can  be 
generated.  The  approximate  Plateau  distribution  Nd  = constant  first  pointed 
rut  by  Caledonia  and  Center  is  seen  to  be  an  approximate  solution  of  Eq.  (6.1) 
with  C > 0. 

Because  A''"''(i)  increases  rapidly  with  i,  the  right-hand  side  of  Eq.  (5.13) 

d 

is  not  uniformly  negligible,  and  the  integration  constant  C is  in  fact  not  a 
constant.  Before  attacking  Eq.  (5.13)  directly,  we  shall  first  introduce  the 
nominal  plateau  distribution  N(i)  defined  as  follows: 


-b5  ^ , 

N(0)  = , 

— 8 9 

where  is  a constant  to  be  determined  later.  Integrating,  we  obtain  ’ 

★★ 

£N(i)  = (i) 


(6.3) 


(6.4) 


where 


.VT  /*i 

N**<n 


(6.5) 


— VT 

and  is  seen  to  be  solely  a function  of  the  VV  and  VT  rate  models.  For  N^>>  c , 

^ wx 

IV(i)  is  approximately  tl  until  the  exponential  ly  increasing  ( i ) intervenes, 

0 a 

bending  it  sharply  downward.  The  knee  of  the  bend  occurs  approximately  at  the 

★★  irk  ^ 

intersection  of  N(i)  and  N (i).  The  location  i where  N'U)  vanishes  is  re- 


lated implicitly  to  by 


^ r 

2b6  J 


A^^di  . 
a 


(6.6) 


For  i > i**,  N(i)  is  identically  zero. 
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We  shall  presently  show  that  N(i)  is  an  approximation  to  N(i)  in  the 
Plateau  region.  Subtracting  Eq.  (6.3)  from  Eq.  (5.13)  we  obtain 


di2 


^nN-6 


1 « 1 < 1 


b6  ^ (N^)^-  Q + e'^'^'g, 


(6.7) 


where 


g(i)  5 (N£  - N£)aJ^ 

a 


(6.8) 


Integrating  Eq.  (6.7)  from  0 to  i**  and  noting  that  the  net  VV  contribution 
must  vanish  identically,  we  obtain 

. ** 

/’  9di- 


‘■“o  = '’total 


-VT 


(6.9) 


which  is  simply  the  quanta  conservation  equation  (3.14;  expressed  in  an  alter- 
native form.  Integrating  Eq.  (6.7)  from  i to  i**,  we  obtain 

2 


where 


b(Nf) 


(N^)' 


^ ^uN-6 


— 2 

b6(!'Itl)  , i « i < i 


(N£) 


b6 


** 


/’«■ 


(6.10) 


(6.11) 


Eq.  (6.10)  is  simply  Eq.  (6.1)  with  a variable  C = b6(N^)^. 


VT. 


We  shall  choose  N(i)  by  the  requirement  that  the  integrals  of  s g(i)  in 
Eqs.  (6.9)  and  (6.11)  be  uniformly  small.  Under  strongly-pumped  conditions, 

Eq.  (6.9)  yields  = ^^total'^^®^^  yields  N(i)  in  terms  of  N(i) 

and  Q(i).  Since  g(i)  contains  Ay"^(i)  which  generally  increases  exponentially 

a 

with  i,  we  conclude  further  that  N(i)  must  approach  fT(i)  rapidly  as  i increases. 
Only  in  the  vicinity  of  i**  where  the  local  values  are  both  small  can  N(i)  begin 
to  deviate  from  N(i)  and  to  approach  a local  equilibrium  Boltzmann  there.  In 
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Other  words,  N(i)  is  the  leading  approximation  to  'i(i)  in  the  Plateau  reoion, 

if  if  irit 

N (i)  is  the  locus  of  the  Plateau-Bol tzmann  knee  of  the  solutions,  and  i 
is  the  location  of  the  Plateau-Bol tzmann  boundary. 

Using  Eq.  (4.1)  for  in  Eq.  (6.9)  and  comparing  with  Eq.  (3.14),  we 

obtain 


= n . n 

m 0 a a 


. ** 

P gd 
Jn 


(6.12) 


The  first  term  represents  the  rate  of  upward  transport  of  quanta  by  nearly- 
resonant  VV  collisions  and  is  seen  to  be  independent  of  VT  rates  even  though 
these  quanta  are  eventually  deactivated  at  high  levels  by  VT  collisions.  The 
second  term  is  a small  correction  representing  the  rate  of  local  low-level  VT 
deactivations  and  can  be  computed  after  M(i)  has  been  found. 


VII.  THE  TREAflOR-PLATEAU  TRANSITION 

★ 

In  the  vicinity  of  the  Treanor-Plateau  boundary  i , the  solution  is  expected 

to  approach  N ('’’^min^  f'(il'iQ)  on  the  other.  The 

value  of  i*  or  i . can  be  related  to  by  a local  transition  analysis  as 
min  0 

f 0 1 1 ows . 

We  shall  tentatively  assume  that  i » i so  that  Eq.  (6.10)  remains  valid 
in  this  transition  region.  Since  the  dominant  external  excitation  is  assumed  to 
occur  mainly  in  the  first  few  levels,  the  Q(i)  contribution  in  Eq.  (6.11)  can  be 
neglected,  yielding  I'J(i)  = N(i)  for  i ^ i*.  For  i < i*,  accurate  representation 
of  N(i)  is  not  required  provided  that  (ri/:l)^  becomes  locally  sm,all.  Thus,  in  the 
transition  region,  N(i)  may  be  approximated  by  a constant,  e N(i  ).  The  equa- 
tion governing  the  Treanor-Plateau  transition  is  then^ 


1 d" 
6 di^ 


TnN  = 1 - 


2 


(7.1) 
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» 


The  magnitude  of  N*  is  assumed  to  be  small  in  comparison  to  the  ground  level 
population 

(N*/Nq)'  « 1 • (7.?) 

= % — 

Even  though  both  fj  = N*  and  the  differential  approximation  underlying  Eq.  (7.1) 

cannot  be  formally  justified  for  0 _<  i = i , it  is  nevertheless  clear  that  the 

TRR 

transition  solution  will  automatically  approach  N (''l^min^  as  i 0 so  long 
as  Eq.  (7.2)  is  satisfied.  Consequently,  we  shall  use  Eq.  (7.1)  in  the  extended 
region  0 £ i = i*  to  describe  the  Treanor-Plateau  transition.  The  values  of  IJ* 
and  i^^^  are  yet  unknown  but  will  be  determined  presently. 

Eq.  (7.1)  can  be  integrated  once  to  yield 


where  N(i  » i*)  = N*  has  been  imposed.  Integrating  again,  we  obtain  an  impMcit 
representation  of  the  Treanor-Plateau  transition  solution:^ 


kinetic  rate  models. 


As  mentioned  previously,  this  transition  solution  automatically  approaches 

TDD  2 

N (i)  as  i 0 provided  (U*/N^)  <<  1.  The  value  of  i^^.^  appearing  in  Eq.  (6.2) 

TRR 

can  now  be  determined  by  matching  the  derivative  of  CkN  and  Cull  at  i = 0,  yield- 
ing 

i 
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,TRR,. 


Since  is  the  location  of  the  Treanor  minimum  of  the  particular 


which  osculates  N(i)  at  i = 0,  we  may  conveniently  identify  it  to  be  the 


Treanor-Plateau  boundary  i*.  Hence,  the  locus  of  the  Treanor-Plateau  boundary 


N*(i)  can  be  obtained  by  inverting  Eq.  (7.5)  and  setting  i^^^  = i*.  V.'e  obtain 


N*(i)  = 0.61  Nq  exp  ( -6iV2 


(7.6) 


For  given  value  of  N^,  the  intersection  of  N (i)  and  lUi;!!^)  immediately  yields 


the  values  of  N*(N  ) and  i*(H  ). 


The  above  analysis  may  be  refined  by  using  a more  accurate  representation 
of  N’(i).  In  addition,  in  the  vicinity  of  i**  a separata  Plateau-5ol  tzmann  trans- 
ition analysis  is  required.  Since  these  refinements  are  not  of  major  practical 
interest,  they  will  not  be  presented  here. 

The  assumption  i >>  i invoked  earlier  can  now  be  examined.  Using  Eq.  (7.5) 


for  i*,  we  obtain  a lov/er  bound  for  the  value  o~ 


VV 

0 » 


1 


(7.7a; 


Furthermore,  in  the  derivation  of  Eq.  (5.6)  which  eventually  led  to  the  differ- 


-VV 


ential  approximation,  contributions  to  F^.  by  non-resonant  VV  collisions,  in- 


cluding those  with  partners  from  the  highly  populated  lower  levels,  have  been 


\J\t 

neglected.  A rough  estimate  yields  the  following  condition  on  A ‘(;;T): 


>VV 


/ ^ -r\  . , bfN*(  i ) . '7^  .* 

(C  = i/l  ;T)  « — a 1 = 1, 


(7.7b) 


which  must  be  satisfied  in  addition  to  Eq.  (7.7a).  If  Eqs.  (7.7)  are  not 
satisfied,  the  Treanor-Plateau  transition  will  depend  on  detailed  information 


VV 

on  A (^');T)  and  must  be  studied  using  the  original  discrete  form  of  the  master 


kinetic  equations.  For  CO,  Eqs.  (7.7)  are  comfortably  satisfied  by  r:ost  VV  rate 
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models  assumed  in  theoretical  calculations. 


VIII.  SUMMARY  OF  RESULTS 

By  taking  advantage  of  the  dominance  of  the  nearly-resonant  VV  collisions, 

we  have  derived  using  the  differential  approximation^ the  loci  of  the  Treanor- 

★ **  * 

Plateau  boundary  N (i)  and  the  Plateau-Bol tzmann  knee  N (i).  N (i)  is  shov/n 

to  be  independent  of  kinetic  rate  information  provided  that  Eqs.  (7.7)  is  satis- 
fied.  N (i),  on  the  other  hand,  is  found  to  depend  on  both  VV  and  VT  rates.  ’ 

At  any  given  temperature,  these  loci  can  readily  be  computed  and  plotted  as  illus- 
trated in  Fig.  1.  The  family  of  nominal  plateau  distribution  I^iill^)  can  also 
be  plotted  and  its  intersection  with  N (i)  yields  immediately  the  values  of 

and  In  the  first  few  levels,  fi(i)  is  osculated  by  a Treanor  dis- 

tribution  iJ  (i,i  . ) with  i . --  i (h').  As  i increase  , N(i)  cresses  N (i) 

mm  mm  o 

via  the  Treanor-Plateau  transition  solution  to  approach  fT(i;fiQ)  in  the  Plateau 

region,  bending  sharply  downward  as  it  crosses  H (i)  and  tending  toward  a local 

★ ★ 

Boltzmann  distribution  near  i (‘1q)-  The  strongly-pumped  condition  is  satisfied 

★ 

whenever  the  solution  passes  substantially  above  the  intersection  of  N (i)  and 

•kic 

N (i)  so  that  a distinct  Plateau  region  exists.  Fig.  1 also  shows  a comparison 

of  the  present  results  with  exact  numerical  solutions  of  the  original  discrete 

equations  for  a CO-He  1:9  mixture  at  100°I(  using  a computer  code  developed  by 
1 5 

Boulnois.  The  agreement  is  surprisingly  good  in  view  of  the  approximations 
invoked.  Note  that  these  discrete  solutions  were  computed  using  39  levels  which 
appeared  to  be  inadequate  for  the  tv/o  more  strongly-pumped  cases. 

The  equivalent  ground  level  vibrational  temperature  T defined  by  fitting 

B ^ 

the  first  few  levels  by  N (T  ) can  readily  be  computed: 


^^ib 

2ar2  ^ - 1 


(8.1) 
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I 


■« 


Other  gross  quantities  of  interest,  such  as  c(N^,T)  and  >’'’)  > can  also 

be  computed.  The  value  of  or  its  time  evolution  is  governed  by  the  quanta 

,VT 


★ ★ 


conservation  equation  (3.14)  with  D given  by  Eq.  (6.12) 

d ^/t7  ri  _ VV  1 ^ gj. 


a(N^,T)  = V f1  - r\J 
dt  0 ' ^ e m 

i=0 


(8.2) 


.VT 


The  small  correction  term  eg  representing  low-level  local  deactivations  may 
be  included  with  no  difficulty.  Once  N^(t)  has  been  found,  the  vibrational 
state  of  the  molecule  of  interest  is  completely  determined.  Hence  under  quasi- 
equilibrium conditions  can  be  considered  simply  as  a state  variable  which 
together  with  T completely  specifies  N^(N^,T),  o(N^,T),  and  T^  (!J^,T).  The 
rate  of  dissipation  of  vibrational  energy  however,  is  not  a state  function 

and  must  be  computed  from  Eq.  (3.17)  instead  of  Eq.  (3.19): 

00 

d 

r.  zit.  - 
i=0 


%ib  = I 1-^i  -|t%ib('^o’T)  • 


(8.3) 


IX.  DISCUSSION 

VV  VT 

The  kinetic  rates  have  been  expressed  in  Section  II  in  terms  of  w , oc,  , 

a 

» a'^'^(5;T),  and  Ay'*'^(T).  The  same  I-  is  used  for  both  VV  and  VT  rates  because 
of  its  theoretical  significance.  Fo"  anharmonic  oscillators,  E-  is  often  repre- 
sented by 


I = 


i+1 

1 -C^ai 


(9.1) 


where  C.|  is  a positive  constant  of  order  unity.  The  VV  adiabaticity  factor 
VV, 

A (^',T)  is  often  assumed  to  be  the  sum  of  two  exponentials  representing  contri- 
butions from  long-range  and  short-range  forces.  As  pointed  out  by  Brau,^  so  long 
as  Eqs.  (7.7)  are  satisfied  the  detailed  shape  of  A ‘(^;T)  is  unimportant.  Only 


1 
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b(T),  its  second  moment,  appears  in  the  theory.  On  the  other  hand,  the  present 

theory  shows  that  accurate  information^^  on  (T),  particularly  at  high  levels, 

I ,a 

VT 

is  required.  Since  A.  ,(T)  is  expected  to  depend  on  i exponentially,  we  may 

1 ,a 

assume  (see  Eq.  (2.9)  ) 


A*^,(T)  = exp 
VT 

where  9 may  be  temperature  dependent.  Consequently,  all 

VT/  > 

mation  can  be  represented  by  only  five  parameters,  uj  (T), 

d 


(9.2) 

the  kinetic  rate  infor- 
uj''''(T),  C^,  b(T),  and 


Using  Eqs.  (9.1)  and  (9.2),  N(i)  can  be  rewritten  as 

i 


N(i) 


1 -C-iCti 
1+1 


^0  ■ ^0  Z 1^2 

j=0  ^ 


(9.3) 


where 


n uj  4ba0^ -L. 

m vib 


, C2(T)  = a6^ 


ib  ■ T • 


(9.4) 


It  is  seen  that  all  kinetic  rate  as  well  as  energy  level  spacing  information 

are  represented  here  by  only  three  parameters,  C^,  C.j , and  C^.  The  simplicity 

cf  Eq.  (9.3)  suggests  the  possibility  of  deducing  the  values  of  C^(T),  , and 

C2(T)  from  a family  of  experimentally  measured  distributions.  At  the  present 

time,  experimental  rate  information  at  high  levels  is  rather  limited;  most  re- 
18  19  20 

ported  data  ’ ’ are  concerned  primarily  at  the  lower  levels.  The  values  of 

Co(T),  C.J , and  C2(T)  deduced  as  suggested  will  be  welcome  complements  to  con- 
ventional measurements. 

Brau^  was  first  to  use  the  differential  approximation,  deriving  Eq.  (5.13) 
and  obtaining  a family  of  approximate  solutions  N(i;C)  from  Eq.  (6.1),  including 
the  Treanor-Plateau  transition  solutions.  The  cumulative  effects  of  VT 
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deactivations  was  neglected  and  the  relationship  between  C and  the  strength  of 

8 9 

external  excitation  was  not  given.  Gordiets  et  al . ’ first  recognized 

the  importance  of  VT  (and  radiative)  deactivations  in  the  Plateau  region,  and 
obtained  approximate  Plateau  distributions  from  Eq.  (5.13)  by  replacing  d^^nN/di^ 
with  -25,  an  empirical  factor  suggested  by  curve-fitting  of  direct  numerical 
solutions.  In  the  present  work,  this  empirical  adjustment  which  is  difficult  to 
justify  mathematically  is  avoided  by  the  artifice  of  introducing  N(i)  by  defini- 
tion. In  their  effort  to  relate  (or  C)  with  > Gordiets  et  al . evaluated 

D from  Eq.  (3.16)  by  omitting  N (i)  from  N(i)  in  the  Plateau  region.  By  noting 


that 


. ★★ 
1 


VT 


I 

j=0 


= 2(b5N 


j=0  j=0 

WT 

it  is  relatively  straightforward  to  show  that  D (see  Eq.  (10)  of  Ref.  8)  is 
approximately  twice  the  correct  result  given  here.  An  added  advantage  of  the 
present  procedure  is  the  simplicity  of  the  expression, Eq.  (6.12),  obtained  for 

showing  clearly  the  contributions  of  high-level  and  low-level  deactivation 

★ ★★ 

contributions.  Together  with  the  explicit  identification  of  N (i)  and  N (i) 
as  the  appropriate  loci,  the  present  theory  yields  a complete  and  quantitative 
picture  of  vibrational  relaxation  problems  under  strongly  pumped  conditions. 

ACKNOWLEDGMENT 

This  research  was  supported  by  AFOSR  (F44620-73-0059) . 


123 


REFERENCES 

1.  A.  I.  Osipov,  Soviet  Phys.  Doklady  102  (1960). 

2.  C.  F.  Treanor,  J.  W.  Rich,  and  R.  6.  Rehm,  J.  Chetn.  Phys.  48,  1798  (1968). 

3.  G.  E.  Calendonia  and  R.  Center,  J.  Chem.  Phys.  ^ 552  (1971  ). 

4.  B.  Abraham  and  E.  R.  Risher,  J.  Appl . Phys.  ^ 4621  (1972). 

5.  J.  W.  Rich,  J.  Appl.  Phys.  2719  (1971). 

6.  W.  R.  Lacina,  Northrop  Corp.  Report  71-32R,  pt.  II,  (1971). 

7.  C.  A.  Brau,  Physica,  ^ , 533  (1972). 


8. 

6. 

F.  Gordiets,Sh. 

S. 

Mamedov , 

and 

L.  A.  Shelepin, 

Ah.  Eksp, 

. Teor.  Fiz 

61. 

, 1287  (1974). 

9. 

B. 

F.  Gordiets  and 

Sh, 

. S.  Mamedov, 

Prik.  Mat.  Teor, 

. Fiz.  3, 

13  (1974). 

10. 

R. 

N.  Schwartz,  Z. 

I. 

Slawsky , 

and 

K.  F.  Herzfeld, 

J.  Chem. 

Phys. 

(1952). 

11.  R.  N.  Schwartz  and  K.  F.  Herzfeld,  J.  Chem.  Phys.  2^,  767  (1954). 

12.  J.  Keck  and  G.  Carrier,  J.  Chem.  Phys.  ^ 2284  (1965). 

13.  L.  D.  Landau  and  E.  Teller,  Physick.  Z.  Sowjetunion  1_0,  34  (1936). 

14.  K.  F.  Herzfeld  and  T.  A.  Litovitz,  Ahscmtion  and  Disperson  of  Ultrascnio 
Waves,  Academic  Press,  N.  Y.  (1959). 

15.  J.  L.  Boulnois,  ’’Carbon  Monoxide  Vibrational  Kinetics  in  a Transverse 
Electric  Discharge  Wave  Guide  Laser",  Ph.D.  Thesis,  Department  of  Aerospace 
and  Mechanical  Sciences,  Princeton  University,  Princeton,  N.  J.  (1975). 

16.  M.  Verter  and  H.  Rabitz,  J.  Chem.  Phys.  M,  2939  (1976). 

17.  W.  S.  Drozdoski,  R.  M.  Young,  R.  D.  Bates,  Jr.,  and  J.  K.  Hannock,  J.  Chem. 
Phys.  65,  1542  (1976). 

18.  R.  C.  Millikan  and  D.  R.  White,  J.  Chem.  Phys.  98  (1963). 

19.  D.  J.  Miller  and  R.  C.  Millikan,  J.  Chem.  Phys.  5^,  3384  (1970). 


124 


FIGURE  CAPTION 


Fig.  1 


Comparison  with  nutrerical  solutions.  The  solid  dots  are 
numerical  solutions  for  a 1:9  CO-He  mixture  at  10C°K 
with  various  pumpings  at  the  ground  level.  The  YV  rates 
are  taken  from  Rich^  and  the  VT  rates  are  recent  theoretical 
data  by  Verter  and  Rabitz.^^  Thirty-nine  levels  are  used 
which  appear  insufficient  for  the  two  more  strongly-pumped 
cases.  The  analytical  solutions  are  shown  as.  solid  lines. 
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